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We analyze two representative systems containing a three-phase-contact line: a liquid lens at a 
fluid-fluid interface and a liquid drop in contact with a gas phase residing on a solid substrate. In 
addition we study a system containing a planar liquid-gas interface in contact with a solid substrate. 
We discuss to which extent the decomposition of the grand canonical free energy of such systems into 
volume, surface, and line contributions is unique in spite of the freedom one has in positioning the 
Gibbs dividing interfaces. Before being able to define the line contribution to the grand canonical 
free energy per unit length, the so-called line tension, first certain properties have to be attributed to 
the bulk phases and to the interfaces; these are chosen in agreement with what is known for spherical 
drops and for planar interfaces. In the case of a lens it is found that the line tension is independent 
of arbitrary choices of the Gibbs dividing interfaces. In the case of a drop, however, one arrives at 
two different possible definitions of the line tension. One of them corresponds seamlessly to that 
applicable to the lens. The line tension defined this way turns out to be independent of choices of the 
Gibbs dividing interfaces. In the case of the second definition, however, the line tension does depend 
on the choice of the Gibbs dividing interfaces. We also provide equations for the equilibrium contact 
angles which are form-invariant with respect to notional shifts of dividing interfaces which change the 
description of the system but leave the density configurations unchanged. In order to reconcile the 
transformation laws for the line tensions under notional changes of dividing interfaces, which follow 
from the decomposition of the grand canonical free energy, with the principle of form invariance and 
moreover with the invariance of observables under notional changes, additional stiffness constants 
attributed to the line must be introduced. The choice of the dividing interfaces influences the 
actual values of the stiffness constants. We show how these constants transform as a function of the 
relative displacements of the dividing interfaces. Our formulation provides a clearly defined scheme 
to determine line properties from measured dependences of the contact angles on lens or drop 
volumes. This scheme implies relations different from the modified Neumann or Young equations, 
which currently are the basis for extracting line tensions from experimental data. These relations 
show that the experiments do not render the line tension alone but a combination of the line tension, 
the Tolman length, and the stiffness constants of the line. In contrast to previous approaches our 
scheme works consistently for any choice of the dividing interfaces. It further allows us to compare 
results obtained by different experimental or theoretical methods, based on different conventions of 
choosing the dividing interfaces. 

PACS numbers: 68.05.-n, 68.08.-p 



I. INTRODUCTION 

The handling of small liquid droplets in contact with 
a gas phase on top of a solid substrate, or of liquid lenses 
at the interface of two other fluid phases plays an impor- 
tant role in the context of microfluidics (see, e.g., Refs. 
[3, 0, 0, Q)- For droplets or lenses with linear exten- 
sions in the nanometer regime a proper thermodynamic 
description requires to account for not only the bulk and 
surface properties but also for the special properties of 
matter in the vicinity of three-phase contact. In order to 
capture the contribution of the three-phase-contact re- 
gion to the relevant thermodynamic potential, a certain 



excess contribution to the appropriate free energy that 
scales with the linear extension of the system is associ- 
ated with the three-phase-contact line, defined by a com- 
mon intersection of the interfaces meeting in the region 
of three-phase contact. 

The liquid lens is an example of an inhomogeneous sys- 
tem in which three thermodynamic phases, say a, (3, and 
7, coexist in a (constrained) equilibrium. The thermody- 
namic coexistence of bulk phases is provided by specific 
choices of the thermodynamic state of the system whereas 
their spatial coexistence follows from appropriate bound- 
ary conditions. We consider the following set-up. A drop 
of the non- wetting (3 phase is placed at the interface be- 
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tween two other fluids. A microscopically thin equilib- 
rium wetting film of (3 phase is formed at the interface 
and for suitably chosen substances and conditions a sur- 
plus of f3 phase forms a lens at the a-7 interface. (In 
principle the exploration of the configuration space al- 
lows for shifting the lens laterally along the interface; 
however, for the purposes of the present paper one may 
disregard this degree of freedom.) For such a system two 
basic scenarios are possible. The lens {(3 phase) can ex- 
change matter with the surrounding phases, so that the 
chemical potentials in all phases are equal. There are 
cases in which the lens is in a stable equilibrium with the 
surrounding phases; in other cases the lens is unstable 
but could be stabilized by imposing suitable constraints. 
Altcrnativly, one may consider a nonvolatile liquid {(3 
phase), i.e., one constrains the volume of the liquid while 
chemical equilibrium is not attained. 
In addition to the concept of interfaces separating the 
bulk phases and of interfacial or surface tensions, the 
systematic thermodynamic description of such systems 
leads to considering the contact line L a p 1 along which 
the three interfaces meet. The line tension t is attributed 
to the contact line C a ^. It is defined as the 'line con- 
tribution' to the grand canonical free energy fl of the in- 
homogeneous system per unit length of the contact line 
S S 0, B Q which is the leading term left after subtract- 
ing volume and surface or interface contributions from Q, 
i.e., 

fi — ( Vlt - A a (3 <J a f3 — A aj a aj — Aflry (Jf3j 



S.l.t. 



(1.1) 



With s.l.t. we denote subleading terms which vanish for 
L a f3-f — ► 00. The symbol ui K denotes the grand canon- 
ical free energy density of the homogeneous bulk phase 
k (k = a,/3,7), V K is the volume assigned to this phase, 
A KK i is the area of the k-k' interface, cr KK < is the corre- 
sponding interfacial tension, and L a fj^ is the length of 
the three-phase-contact line. The definition of t in Eq. 
(jl.ip refers to a reference state in which uniform bulk 
phases are extrapolated right up to the interfaces, and 
analogously the interface or surface properties of later- 
ally homogeneous interfaces or surfaces are extrapolated 
right up to the contact line £ Q /3 7 . In Eq. (jl.lj) we do 
not take into account contributions to f2 related to the 
presence of walls enclosing the whole system; the only 
inhomogeneities of the system which are relevant in the 
present analysis are those related to the spatial coexis- 
tence of the three thermodynamic phases a, f3, and 7. 
The above decomposition allows one to calculate the line 



tension for a given thermodynamic system provided the 
quantities on the lhs of Eq. (jl.ip are determined in sep- 
arate preceding steps involving similar considerations of 
suitable thermodynamic limits. 

A different typ of an inhomogeneous fluid system is 
a sessile drop on a solid substrate. In this case only 
two thermodynamic phases coexist whereas the substrate 
acts as an inert spectator phase. In such cases the three- 
phase-contact line corresponds to the region where the 
interface between two coexis ting thermodynamic phases 
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In the literature the term 
line energy or line elasticity is also frequently used in 
connection with an extra surface energy associated with 
the deformation of the three-phase-contact line as the 
result, e.g., o f interactions with surface defects (see, e.g., 
Refs. (TflElIllIIEIiHl)- This quantity has to be clearly 
distinguished from the line tension as introduced via Eq. 

In the literature the notion of line tension is also used 
to describe the one-dimensional interface of two coexist- 
ing, intrinsically two-dimensional phases such as liquid- 
and vapor-like phases in Langmuir-Blodget films. This, 
however, corresponds to the lower-dimensional version of 
an interfacial tension and not to the coexistence of three 
phases as considered here. 

There is a growing body of literature describing both ex- 
perimental and theoretical investigations of the line ten- 
sion and already a number of reviews have been published 
on that subject 0, E3, H3, [13, E3] • Experimental in- 



vestigations were carried out for drops on solid substrates 
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511 ], for liq uid lenses at the interface 



of two fluids 

spherical particles or bubbles at the interface of two flu 
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5J, |55j, [5J, |57J, and for 
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line tension in epitaxial growth has been discussed [68 
On the theoretical side the theory of capillarity has been 
extended in a phenomenological way by taking into ac- 
count line contributions and the consequences of these 
extensions have been explored [ S t 1 18l l36l . I52L |6(1 \69L 170 . 
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Furthermore there are microscopic calculations of the line 
3, 0, 0000 [93, 



tension 




96J and many studies, in most cases based on mi- 
croscopic theories as well, concentrate on the behavior 
of the line tension in the v i cinity of a we t ting transition 
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Theoretical studies are also devoted to the examination 
of th e influence of the droplet size on the line tension 
11081 and of the electrostatic contributions to the line ten- 
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sion which arise in the presence of surface charges 109 1 . 
The boundary tension between two coexisting wetting 
films of di fferent th i ckne sses was studied as well (see, e.g., 
Refs. 93Ul06l . lllOUlll| ). A drop-size dependence of the 



contact angle was studied in molecular dynamics simu- 
lations wit hout deducing values for the line tension (see, 



e.g., Ref. |ll2| ). Finally, the line tension was studied 
also via m olecular dynamics simulations (see, e.g., Refs. 



113l 114 1) or via an analysis of probability distributions 



as obtai ned f rom Monte Carlo simulations at three-phase 
contact 
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A closer inspection of these results reveals that certain 
aspects of the line tension either give rise to conflicting 
statements or remain unaddressed, leading us to set out 
to clarify the following basic questions: 

• In which sense is the concept of the line tension 
well defined? 

• Is it sufficient for the determination of equilibrium 
shapes to characterize the thermodynamics of the 
contact line by a line tension only? 

These questions arise because it is not obvious that the 
definition of r via Eq. (jl.ip leads to a unique result. 
The reason lies in the arbitrariness in the definition of 
the position of the interfaces between adjacent phases. 
This arbitrariness is due to the smooth spatial transition 
between the adjacent phases. Once the density distribu- 
tions of the fluid components across the interfacial region 
are known, for example on the basis of scattering experi- 
ments or atomic force microscopy 4CJ, |46| . |47| , or theoret- 
ically by simulations or density-functional calculations, a 
criterion has to be applied which fixes the interface posi- 
tion somewhere in the transition region. However, there 
exists a multitude of sensible choices. The arbitrariness 
in the definition of the interface positions leads to an arbi- 
trariness in the definitions of the volumes, areas, lengths, 
and to s ome extent even interfacial tensions (see, e.g., 
Refs. @,[Ill[llJJ^[l3,[l2B,[l2ll). Although it was 
shown by Widom 0,[8 1 that r is independent of the par- 
ticular choice of the dividing interfaces in the case of a 
straight three-phase-contact line in a system with three 
thermodynamically coexisting fluid phases separated by 
planar interfaces, it is not clear that the same will be true 
in other cases as well, e.g., if two fluid phases are in con- 
tact with an inert solid phase, or for systems in which the 
curvatures of the interfaces and of the contact line plays a 
role. A first discussion, which naturally raises the issues, 
concerning the uniqueness of the line tension for the case 
in which two fluid phases are in contact with and meet 
at an inert solid phase, is given in Sect. 3. This discus- 



sion, however, turns out to be incomplete and it leads to 
a contradiction with one of the results obtained in Sect. 
5, which states that the line tension as defined there is 
unique. A thorough discussion and the resolution of the 
contradiction is given in Sect. 6. In Sect. 5 we mainly 
investigate the three-phase-contact line in systems with 
curved interfaces. On the other hand, just these systems 
are studied experimentally if one attempts to determine 
the value of the line tension, e.g., via measuring the de- 
pendence of the contact an gle on droplet size (see, e.g., 

Refs. [ii SSH, 0, HHliiliilii.liEy,^). 

The curvature dependence of interfacial tensions leads 
to contributions to the free energy which scale with the 
linear extension of the system, i.e., very much like the 
line contribution. Since in many cases interface curva- 
ture and the length of the three-phase-contact line cannot 
be varied independently, these two contributions are not 
distinguishable per se. Thus, line tension and curvature 
effects on the interfacial tensions have to be discussed si- 
multaneously. Moreover, the curvature expansion of the 
surface tension has to be known in advance before the 
line tension may be determined from Eq. (|1.1|) . 
Since in the present context it is unavoidable to con- 
sider curvature effects on the interfacial tensions and be- 
cause it is known that a consistent description of curved 
(spherical) interfaces must take into account — for gen- 
eral dividing interfaces — derivatives of the surface ten- 
sions with respect to their radii of c urvature, i.e., their 
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1341 I135L Il36j ). we expect that a com- 
plete thermodynamic description of three-phase-contact 
lines requires the introduction of further material param- 
eters in addition to the line tension. Indeed, a num- 
ber of theoretical analyses have appeared in the litera- 
ture in which additional material parameters were intro- 
duced in order to characterize the three-phase-contact 
line 0, Izl, S S Izi- Re cent studies have been car- 
ried out by Rusanov et al. |l37| for a drop on a solid 
substrate, and by Widom and coworkers in connection 
with the line analog ue o f the Gibbs adsorption equation 
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143j for a straight contact line 
at a genuine three-phase contact. In Sect. 5 we sys- 
tematically study these issues guided by the concept of 
form-invariance of the basic equations under so-called no- 
tional shifts of dividing interfaces (for a definition see, 
c.f., Sects. 2 and 5). 

As a prerequisite to a thorough study of these issues 
we also have to answer similar questions related to the 
concept of surface tensions. However, before discussing 
them and in particular the main questions concerning 
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the line tension, we point at related issues recurring in 
the literature which can be judged only after the basic 
questions raised above have found a satisfactory answer. 
These issues are: 

• What is the typical magnitude of the line tension ? 
Of course, the line tension depends on the thermo- 
dynamic state of the system and the number of rele- 
vant thermodynamic degrees of freedom varies from 
system to system. For example, for certain systems 
the temperature dependence of the line tension be- 



sitions la . 



17, 


18 




20 


21 


10( 


. 


101 




102 



22, 33 



103 



9(1 



91 



92, 93, 



On the other hand, for many systems away from 
such special thermodynamic states there is now 
widespread agreement that the value of the line 
tension is of the order of 10" 11 N 0, 0, 0, Efi , 
13]. With the experimental accuracy available at 
present the corresponding prefactor distinguishes 
between different systems. 

• What is the sign of the line tension? 
Various experimental and theoretical investigations 
lead to values of the line tension which incl ude both 
signs EH Ezl, El S y , EH, EI3- Since - 
contrary to the interfacial tension - there exists no 
thermodynamic argument for a specific sign of the 
line tension 0, 148 1. so that experimental findings 
for the line tension with a certain sign cannot be 
discarded from the outset. (Even in this respect a 
cont rary statement may be found in the literature 
|l49| .) On the other hand, if one is interested in 
the temperature dependence of the line tension for 
a specific system which undergoes a first-order wet- 
ting transition, both theoretical and experimental 
results show that upon increasing the temperature 
towards the wetting transition temperature the line 
tension change s sig n from negative to positive val- 



ues 



17|, [20j, l33|, [33, [44 1 . 



• Is a negative line tension compatible with the struc- 
tural stability of drops? 

The mesoscopic analysis of line tensions including 
small wavelength fluctuations confirms that nega- 
tive values of the line tension do not lead to insta- 
bilities 
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• What does the line tension depend on? 
Similarly to the interfacial tension it is a function 
of the thermodynamic state of the system. For 
straight contact lines, defined by intersecting pla- 
nar interfaces, the phases involved must be at sta- 



ble thermodynamic coexistence. For example, one- 
component fluids in contact with an inert substrate 
must be at liquid-vapor coexistence /i = ^q{T), 
where /i denotes the chemical potential. Thus in 
this case the line tension is a function of tempera- 
ture only. If the fluid in contact with a substrate 
consists of a binary mixture of A and B particles, 
the fluid must be at fluid-fluid coexistence, i.e., 
/i A = /i A (/i B ,T) which leaves two thermodynamic 
variables free. For three-phase contact among three 
fluid phases the system must be at the triple line 
(jj, a = M A (T), ^ B = /i B (T)) of A-rich liquid, B-rich 
liquid, and vapor coexistence so that in this case the 
line tension is again a function of temperature only. 
This dependence can be reparametrizcd in terms 
of the temperature dependent contact angle 0(T). 
The situation is somewhat different for curved con- 
tact lines defined by intersecting curved interfaces. 
Droplets of finite size residing on a substrate or 
liquid lenses formed at planar fluid-fluid interfaces 
are examples in which the contact lines are curved. 
Due to their curvature droplets or lenses remain 
in a (constrained) equilibrium with their surround- 
ing phases, which takes place at chemical poten- 
tials slightly off their values /i l at stable thermo- 
dynamic coexistence with planar interfaces. Under 
these conditions the pressure is different inside and 
outside the droplet or lens. This pressure difference 
and at the same time the size of the droplet or lens 
are prescribed by the chosen chemical potentials or 
alternatively by the temperature and the ambient 
pressure which deviates from that for stable coex- 
istence at the same temperature. Thus in principle 
the line tension may now depend on a further ther- 
modynamic variable in addition to T. This depen- 
dence can be reparametrized in terms of the size of 
the droplet or lensc. 

In addition, from the theoretical point of view the 
line tension r is a functional of both the interac- 
tion potentials between the fluid particles and the 
substrate potential. If the microscopic forces are 
too long-ranged r becomes ill defined while the cor- 
responding interf acial and surface tensions retain 
their validity 17, 150j . i.e., in this case the size de- 
pendence of f2 cannot be described in terms of a 
bulk, surface, and line contribution with a line ten- 
sion which is size- and shape-independent in the 
thermodynamic limit. 

It appears that in extracting line tensions from ex- 
perimental data so far the curvature dependence of 
the interfacial tension, characterized by the Tolman 
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been completely neglected. Therefore the question 
arises to which extent the experimental determina- 
tion of the line tension is affected by the Tolman 
length. 

• Can one measure the line tension? 
In the last chapter we shall point at possible diffi- 
culties to determine line tensions uniquely. 



II. THE SUBSTRATE-FLUID SURFACE 
TENSION 



In the following we discuss a possible source of ambigu- 
ity in the definition of the line tension (see Eq. (jl.lj) ) in 
the case that one of the phases considered in Sec. 1 is an 
inert substrate, i.e., we consider two-phase coexistence in 
the presence of an inert wall rather than coexistence of 
three genuine thermodynamic phases. We focus on the 
issue of non-uniqueness of the liquid-substrate or gas- 
substrate interfacial tensions. To this end we first briefly 
mention the related question of the thermodynamic de- 
termination of the interfacial tension o KK i in a system 
in which two coexisting fluid phases, say k and k' , meet 
along a planar interface 0, @| ■ Using the formula 



lim 



O - V K u K - V K , uj k > 



(2.1) 



the question appears whether the corresponding value of 
the interfacial tension depends on the arbitrary choice 
of the position of the fluid-fluid interface and the corre- 
sponding volumes V K and V K i. In the case of a planar 
interface the value of the interfacial tension does not de- 
pend on the position of the interface, i.e., Eq. (|2.ip leads 
to a unique result. The reason is that the total volume 
V K + V K i as well as the surface area A KK > are independent 
of the location of the dividing surface; in addition due to 
thermal equilibrium one has lu k — u> K >. (For a complete 
discussion of this issue see, e.g., Refs. 0, Q.) 
However, in situations in which one of the phases, say 
phase 7, is an inert substrate with a planar surface, a 
difference to the case of fluid-fluid interfaces arises be- 
cause the inert substrate is not one of the thermody- 
namically coexisting phases. Although it is not a priori 
obvious whether one should treat the substrate as a part 
of the system or not, it is usually considered as an exter- 
nal object which just provides a steep external potential 



defining the boundaries of the system. 
It might appear that the freedom in positioning the di- 
viding interface is less obvious for the solid-fluid interface 
than for the fluid-fluid interface because a solid surface 
is defined rather sharply, say by the positions of the nu- 
clei of the atoms forming the topmost layer. However, 
even then the position at which an actual experimental 
technique (e.g., atomic force microscopy, optical meth- 
ods, etc.) will locate the surface of the solid will certainly 
deviate from the definition given above, because of the 
smooth decay of the substrate potential and of the finite 
extension of the electron cloud of the substrate, and be- 
cause the fluid phase in contact with the solid only grad- 
ually attains its bulk properties. Depending on the kind 
of experiment and on the way the data are analyzed, rel- 
ative shifts in the location of the dividing surface which 
are of the order of one or even several atomic radii are 
conceivable. This shift multiplied by the surface tension 
under consideration yields a force which is comparable 
with the magnitude of the line tension. Moreover, a solid 
substrate in contact with a vapor phase might be cov- 
ered with a thin liquid-like wetting film which is in ther- 
mal equilibrium with the bulk vapor phase. In a ther- 
modynamic description in terms of interfacial and line 
tensions as considered here, this thin wetting film is not 
treated as a separate entity but as part of the solid-vapor 
interface and as such it contributes to the actual solid- 
vapor (gas) surface tension a sg . Accordingly the question 
arises, where to place the solid-vapor interface. It could 
be somewhere in the transition region from the liquid- 
like film to the vapor or at the transition region from the 
solid to the liquid-like film. It is very likely that differ- 
ent experimental techniques imply different conventions. 
(If nonequilibrium situations are considered a sg may be 
different from its equilibrium value if the aforementioned 
thin wetting film - present at thermal equilibrium - has 
not yet formed.) 

In order to explore the consequences of the freedom in 
choosing the dividing surface, we consider a system in 
which a planar inert substrate is exposed to the fluid /, 
i.e., the liquid or gas phase (see Fig. 1). If we consider the 
solid as a part of our system the value of the substrate- 
fluid surface tension a s f follows from the equation 



<7 s f = lim 

Vl, sl A 8 f->-oo 



fi — Vf U)f — V s U>s 



A 



*f 



(2.2) 



where denotes the grand canonical potential of the fluid 
plus that of the substrate, Vf is the value of the fluid vol- 
ume compatible with the chosen location of the planar 
substrate-fluid interface, V s is the corresponding volume 
of the solid, and uif and tu s are the grand canonical free 
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FIG. 1: A fluid phase (f) in the presence of a planar substrate 
(s). Two choices for the position of the planar substrate-fluid 
interface are marked as (1) and (2); they are shifted with 
respect to one another by a distance Sh. In each case the area 
A s f of the substrate-fluid interface is the same. 



energy densities of the fluid and solid, respectively. A s f 
is the area of the planar solid-fluid interface. If the sub- 
strate is in a constrained equilibrium (e.g., no interdiffu- 
sion of solid and fluid particles), one has ujf =/= lu s . 
Since the position of the substrate-fluid interface is not 
unique we can consider another position of this interface 
which is parallel to the first one and located a distance 
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(1) 



-Sh A* 



Sh above it (see Fig.l). In this case V* — V^- —oni\ s f, 

V s {2) = Vs W +ShA sf , wher eas the surface area A s f does 
not change upon the vertical shift of the interface. It fol- 
lows from Eq. (|2.2[) that the values of the substrate-fluid 
surface tensions corresponding to these two choices differ 
by 



r (2) 



(ujf — lo s ) Sh 



(2.3) 



Thus a different choice of the location of the substrate- 
fluid interface corresponds to a redistribution of the free 
energy between the bulk and the surface terms, and to 
different values of the substrate-fluid surface tension. 
Note that in the case of thermodynamically coexisting 
liquid and gas phases the grand canonical free energy 
densities of the two phases (i.e., the negative pressures) 
are equal (u>i = uj g ) and thus the well known conclusion 

7 lg 



a$ follows (see Refs. @, H))- At the same time 



19 ... . 
we emphasize that the difference a sg — o s i (which, e.g., 

enters into Young's law for the contact angle (see, c.f., 
Eq. (|3.2p ) does not depend on this choice of the divid- 
ing surface provided the dividing interfaces between solid 
and gas on one hand and between solid and liquid on the 
other hand are chosen to be at the same height above the 
solid, and provided the solid is in the same state in both 




FIG. 2: Coexisting gas (g) and liquid (1) phases in the pres- 
ence of a planar substrate (s). The planar liquid-gas interface 
meets the substrate with a contact angle 6. Far to the left 
(right) there is the substrate-gas (-liquid) interface. Two 
choices of the position of the substrate-fluid interfaces are 
marked by (1) and (2); their distance is denoted by Sh. This 
results in two parallel contact lines the cross sections of which 
are indicated by the dots. AA\ g and Ay4 s i denote the corre- 
sponding changes in the liquid-gas and substrate-liquid inter- 
facial areas, respectively. For both choices the contact angle 
6 is the same. 



III. THREE-PHASE-CONTACT LINE AT 
INTERSECTING PLANAR INTERFACES: 
UNIQUENESS OF THE LINE TENSION 



In the case of three genuine thermodynamically coex- 
isting phases the issue analogous to the one raised in 
Sect. 2 is whether the value of the line tension t de- 
termined from Eq. Ijl.ip depends on the location of the 
dividing interfaces which in turn determine the location 
of the contact line. One can show [7| that in the case that 
these phases meet along a straight line parallel shifts of 
this line do not influence the value of the line tension. 
In order to study whether the location of the substrate- 
fluid interface affects the value of t, we investigate the 
following configuration of two-phase coexistence in the 
presence of a planar solid substrate. Due to thermal 
equilibrium of the gas and liquid phase one can impose 
lateral boundary conditions such that far to the left the 
substrate is exposed to the gas phase whereas far to the 
right the substrate is exposed to the liquid. This enforces 
the formation of the liquid-gas interface which meets the 
substrate with a contact angle 9 (see Fig. 2). On a 
macroscopic scales the liquid-gas interface is also planar. 

We again consider two parallel positions of the substrate- 
fluid interface at a distance Sh from each other. This re- 



suits in two corresponding, parallel contact lines. From 
simple geometrical considerations it follows that the val- 
ues of the line tensions corresponding to these choices 
differ by 



I (!) W\ 
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(3.1) 



In the calculation leading to Eq. (|3.1|) we have used the 
aforementioned result <Tsg — = trij and the fact 
that the freedom in the choice of the substrate-fluid in- 
terface position does not influence the value of the liquid - 
gas interfacial tension a^ a , i.e., offi = <7 ; ^ = ai g . After 



using Young's equation 



CTsq = Osl + <Jlq COS U 



and identifying 8 — 8q one obtains 



(3.2) 



(a) 



(b) 



a 
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FIG. 3: (a) A liquid lens (phase /3) at the a-7 interface be- 
tween two fluid phases a and 7; r is the radius of the circular 
intersection between the two spherical caps forming the lens. 
The various contact angles are denoted also as a, j3, and 7. 
(b) A sessile liquid drop (phase /3) in contact with its vapor 
(phase a) and a planar substrate (7) forming a contact angle 
9. The drop is a spherical cap. 



IV. LINE TENSION AND CONTACT ANGLES: 
PROBLEMS ASSOCIATED WITH THE 
STANDARD MODIFIED NEUMANN AND 
YOUNG EQUATIONS 



r (2) _ r (l) _ 5h sin6 , Q 



(3.3) 



This result reflects the fact that there is no explicit force 
balance perpendicular to the solid surface. Taken at 
face value the above result would show that the freedom 
in positioning the substrate-fluid interfaces with the 
consequential shift of the contact line is reflected in the 
change of the value of the line tension. This change 
is proportional to the distance Sh between the two 
arbitrarily selected positions of the substrate-fluid inter- 
faces. Numerically, the rhs of Eq. (|3.3|) is comparable 
with t. However, Eq. (|3.3[) is in conflict with a result 
which will be obtained in Sect. V of the present work, 
namely the invariance of the line tension with respect 
to notional changes of the system. This puzzle will be 
resolved in Sect. VI. Here we only note that it turns 
out that the two conflicting statements are based on two 
different definitions of a line tension both of which seem 
to be absolutely compelling from the point of view of 
how they are introduced. A relation between the two 
definitions of t will be given. But we also point out 
that the way how Eq. (|3.3[) has been obtained above 
has to be scrutinized carefully, because, while deriving 
Eq. (|3.3[) . all contributions to SI which are generated by 
separating a subsystem from its surrounding and which 
are proportional to the linear extension of the system, 
have been disregarded. We shall show that in general 
this is not permissible. 



We start this section by considering genuine three- 
phase coexistence in which a lense consisting of a fluid 
phase (3 is located at the interface between two fluid 
phases a and 7 (see Fig. 3(a)). The lense is taken to 
be formed by two spherical caps of different radii inter- 
secting along a circle of radius r. The three-phase-contact 
line of circular shape is accompanied by the line tension r. 
To simplify the notation also the corresponding contact 
angles are denoted by the same symbols as the phases, 
i.e., a, /3, and 7, where a + /3 + 7 = 2-/T. In the absence of 
the line tension the contact angles an, (3q, and 70 fulfill 
the equation (see, e.g., Ref. Q) 



a ai + a a ft cosao + 0737 cos 70 







(4.1) 



If a line-tension contribution to the constrained grand 
canonical free energy SI is included one obtains, from the 
minimization of SI at a constant volume of the liquid 
phase (3, the modified Neumann equation (see, e.g., Refs. 

mQ) 



Cq/3 (cos a — cosao) + 0737 (cos 7 — cos 70) 



r 
r 



(4.2) 

provided neither bending rigidities of the interfaces nor 
further properties attributed to the contact line other 
than the line tension (such as rigidities against changes 
of contact angles) are taken into account. Equation 
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can be equivalently rewritten as 

a a sin A) t 

cosp = cospo — — . (4.3) 

sin ao a a f3 r 

In particular, if one of the phases, say phase 7, is taken 
to represent an inert substrate with planar surface (/3 = 
7T — a ) (see Fig. 3(b)), the above equation turns into the 
modified Young equation (see, e.g., Refs. 32, 71, 72j|) 



cos6> = cos 9 , (4.4) 

where r denotes the radius of the circular substrate- 
phase-a-phase-/3 contact line. Equations (|4.3|) and (|4.4[) 
represent asymptotic formulae valid in the limiting case 
of large lenses or drops, i.e., r governs the leading behav- 
ior for r — > 00. 

We note that in many experiments involving sessile liquid 
drops a so-called line tension r is deduced from measure- 
ments of the contact angle 9 as a function of the radius r 
via fitting the modified Young equation (Eq. (|4.4jl ) to the 
data 0, EH, 37, 40, 46 1. Similar experiments have been 
carried out with lens- like objects 52, 53, 54, 5^, 5^, 57 1. 
However, a closer look at the procedure described above 
of determining r reveals problems in using the modified 
Young (Eq. (|4.4p ) or Neumann equations (Eqs. (|4.2[ 
14. 3j) ) which are related to the freedom in positioning the 
dividing interfaces. A shift of the solid-liquid dividing 
interface (in the case of a sessile drop) by Sh or a change 
of the radius of a spherical interface (for the lens or the 
drop) by SR leads to changes of the contact angles as 
well. From simple geometrical considerations one finds 
that the corresponding changes in cos 9 or cos (3 are of the 
order Sh/r or SR/r, i.e., they are of the same order as the 
corrections stemming from the presence of the line ten- 
sion. Therefore, upon applying the modified Neumann or 
Young equation to the same physical object but with dif- 
ferent choices for the dividing interfaces one would have 
to introduce two different and suitably chosen values of 
r for different dividing interfaces in order to obtain the 
correct relations between the two corresponding contact 
angles. (In the line-tension related correction terms in 
Eqs. (|4.31 14. 4[) the quantities other than r are cither 
independent of the choice of dividing interfaces or their 
changes with the dividing interfaces give rise to higher 
order corrections.) On the other hand, by decomposing 
n in two different ways we had found that for a straight 
contact line at a genuine three-phase contact the line ten- 
sion is independent of the choice of dividing interfaces. 
This contradicts the result of the previous argument. Fi- 
nally, we can also look at what happens if we shift the 
substrate-fluid interface for the sessile drop. If we com- 
pute the difference t^> — enforced by the geometrical 



relations between and 9^ together with the require- 
ment that both 9^ 2 > and O^ 1 ' fulfil the modified Young 
equation we obtain — = — eri g <5hsin0o which has 
the same structure as the value given in Eq. (|3.3|) (ob- 
tained from a comparison of two different ways of decom- 
posing Q), but it has the opposite sign. 
Summarizing these two findings it turns out that the re- 
lations between two line tension values obtained from two 
different decompositions of f2 for two different sets of di- 
viding interfaces are at variance with the currently used 
modified Young and Neumann equations combined with 
elementary geometrical considerations. 
In order to find equations free from the above inconsis- 
tencies we shall investigate in detail two representative 
systems: a liquid lens at the interface between two fluid 
phases and a sessile drop on a substrate. We shall include 
from the outset the effects of curvature on the interfacial 
tensions and we shall explicitly state all conventions used 
in defining different sets of dividing interfaces, and list all 
the properties which are assigned to the interfaces and to 
the reference bulk phases. 



V. A CLOSER LOOK AT LENSES AND DROPS 



In this section we first study a lens-shaped fluid phase 
f3 located at the planar interface between two other fluid 
phases a and 7, and secondly a sessile liquid drop of fi 
phase in contact with its vapor a on top of an unde- 
formablc inert solid substrate 7. 



A. General considerations 

In this context two main questions are addressed. 
First, how does the line tension depend on 'parallel' dis- 
placements of the positions of the dividing interfaces be- 
tween two phases? Secondly, what can be learned from 
the fact that thermodynamic potentials of the total sys- 
tem must be independent of arbitrary conventions re- 
garding the choice of dividing interfaces? In particular, 
we are interested in the resulting requirements with re- 
spect to the structure of the equations relating the con- 
tact angle(s) to the system size, e.g., to the radius r of 
the three-phase-contact line. These questions are moti- 
vated by the inconsistencies encountered for the modified 
Young and Neumann equations in the previous chapter. 
First, we specify the systems in more detail and give the 
rules that we have chosen to define the dividing inter- 
faces. These rules still admit parallel shifts of the divid- 
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Fig. 4: A liquid lens at the planar a-7 interface (horizontal 
lines). Two choices for the interfaces are marked by (1) and 
(2). The relative concentric shifts of the spherical interfaces 
are characterized by [d_Ri] and [dife]- The planar a-7 inter- 
face is shifted together with the spherical ones, as indicated 
by the dashed and solid horizontal lines and the change in the 
amount of 7 phase shown in lighter gray. The corresponding 
contact angles are denoted by a^ 1 ', 7^ 1 - ) , = 2ty — a' 1 - 1 — 7' 1 ' 
and c^ 2 ' , r y^ 2 \ (3^ = — —^ 2%> ; the contact-line radii are 
r' 1 ' and r^ 2 \ Mi and M2 are the centers for the radii of cur- 



vature R\ x >, RY> = R^' + ldRi] and R 2 L > , R [ 2 > = R 2 L, + [dR 2 } 
respectively. 



ing interfaces. We also introduce the quantities that are 
used in order to describe the shifted dividing interfaces. 
It is further necessary to fix the properties attributed 
to the (reference) bulk and the interfaces present in the 
lens or drop, and we do this in accordance with what is 
known about the bulk fluid and the interfacial properties 
of spherical drops. Only after this step is completed the 
properties attributed to the line may be extracted. 
The lens and the drop are completely characterized once 
the density distributions are known for all constituent 
species at given thermodynamic conditions. Only for 
small drops the line tension and curvature effects are ex- 
pected to become relevant, and in this case gravity can be 
neglected. (The validity of Young 's law in the presence of 
gravity is discussed in Ref. 155].) We also assume that 
no other external bulk forces act on the systems. Under 
these conditions, provided the volume of the (3 phase is 
sufficiently large, one can expect that there exist regions 
of the lens (at some distance away from the three-phase 
contact lines) where the density distributions exhibit - to 
a good approximation - radial symmetry relative to one 
of the two centers of curvature Mi and M2; an additional 
'center' at infinity characterizing the planar interface is 
present in the case of a lens (see Fig. 4). In the case 




Fig. 5: A sessile liquid drop on a planar substrate. Two 
choices for the liquid-vapor (/3-a), the substrate-liquid (7- 
P), and substrate-vapor (7-a) interfaces are denoted as (1) 
and (2); the corresponding concentric shifts are characterized 
by [dR] and [dh]. The corresponding contact angles are de- 
noted by and 6^ 2 \ and r^ 1 ' and are the corresponding 
contact-line radii. M is the center of the radii of curvature. 



of the drop we correspondingly expect 'radial' symme- 
tries around one center M plus an additional 'center' at 
infinity characterizing the substrate-fluid interface (see 
Fig. 5). For the lens, the center Mi characterizes the 
density distributions at the a-/3 interface, the center M2 
characterizes those at the /3~7 interface. In view of these 
radial symmetries it makes only sense to consider con- 
centric shifts of the interfaces with respect to the fixed 
centers. The two phases a and 7 are assumed to be sepa- 
rated by a planar interface. If the a-/3 and interfacial 
structures do not overlap except near the three-phase re- 
gion, the interior of the lens is occupied by an almost 
homogeneous (3 phase. This homogeneous (3 phase inter- 
polates smoothly between the radially symmetric density 
distributions around the two centers associated with the 
two interfaces. Similar considerations apply to the drop 
except that the /3~7 interface is planar which can be re- 
garded as a limiting case of a spherical interface with 
infinite radius of curvature. 

From the previous remarks it follows that to a large ex- 
tent the isodensity contours of the a-(3 and (3—j inter- 
faces are segments of spherical surfaces. In order to define 
Gibbs dividing interfaces separating two adjoining phases 
we use the spherical parts of the isodensity contours in 
a two-phase region and extrapolate them into the three- 
phase-contact region where surfaces of constant densities 
actually are no longer spherical. Which of the infinite 
number of surfaces of constant densities is chosen in or- 
der to construct a Gibbs dividing interface is of course a 
matter of convention. 

For the lens, once the centers Mi, Mi are given and the 
radii R\, R2 are chosen according to a certain conven- 
tion we can define a three-phase-contact line of circu- 
lar shape by the intersection of the two spheres {M\,Rx) 
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and (M2,i?2) as indicated in Fig. 4. The third, pla- 
nar Gibbs dividing interface between the phases a and 
7 is placed then in such a way that it coincides with 
the plane determined by the previously defined circular 
three-phase-contact line. Again this choice is a mere con- 
vention but deviating from it would create three different 
lines of intersection between three pairs of interfaces (a- 
f3 intersecting with /3-7, a-(3 with a-7, and a-7 with 
/3— 7). Furthermore, for such a deviating choice a ring- 
shaped volume, with triangular cross section defined by 
these three lines and the connecting interfaces, could not 
be assigned to any of the three phases and would have 
to be treated separately as a new "line phase". This 
would create an unrewarding complication because the 
central idea of introducing the mathematical Gibbs di- 
viding surface is to allow for a natural extrapolation of 
the properties of the adjacent bulk- like phases from both 
sides right to the dividing surface were they are assumed 
to change discontinously. Within this approach there is 
no room for additional phases. A consistent description 
of the actual system follows by defining appropriate sur- 
face excess quantities and their densities, like the surface 
tension. We stick to Gibbs' idea and we use the con- 
vention described above within which all three lines of 
intersections coincide, and no volume filled with a phase 
of an unassigned character is left over. 
The procedure for the drop is analogous. The sphere 
(M, R) defines the a-(3 Gibbs dividing surface. Again, 
the spherical part of the isodensity contour is extrapo- 
lated into the three-phase-contact region. In the same 
spirit as above the a-7 and the Gibbs dividing in- 
terfaces are chosen to lie in the same plane. The three- 
phase-contact line is defined by the intersection of the 
a—f3 Gibbs dividing surface with the common a-7 and 
/3~7 plane (see Fig. 5). Using this convention has the 
advantage of being in agreement with the one chosen for 
the lens in the limit of one of the curvature radii becom- 
ing infinite. Furthermore, this construction creates only 
one line of intersection instead of two. 
Once the dividing interfaces are defined the total volume 
V is subdivided into domains assigned properly to the 
phases a, (3, and 7. The grand canonical potential of 
the system is decomposed into the corresponding bulk, 
intcrfacial, and line contributions. For the lens one has 

Q = - ^2 PkV k + A a pa a p + Ap 1 crp 1 

+ (A - irr 2 ) a al + 2irrr , (5.1) 

where p K is the pressure in the bulk-like phase k (= a, 
(3, or 7) and V K is the corresponding volume. A is the 
area of the planar a-7 interface in the absence of the 



lens. The values of V K depend on the way the dividing 
interfaces are chosen. The total volume V = J2 K ^« °f 
the system is independent of this choice, and independent 
of the physical size of the lens or the drop. 
For the drop one has 

fi = 2J -PnV K + cj 7 U 7 + A af3 a al 3 + TTr 2 ap 1 

k={q,/3} 

+ (A - irr 2 ) er Q7 + 2irrr , (5.2) 

where A a p, Ap~ (= irr 2 for the drop) are the areas of 
the a-(3 and the /3-7 interfaces, respectively. A is the 
area of the planar a-7 interface (in the absence of the 
drop) and it is independent of the choice of the dividing 
interfaces. The area irr 2 denotes the a-7 interfacial area 
replaced by the drop. The radius of the circular three- 
phase-contact line is denoted by r and a a p , ap-y , and a aj 
denote the surface tensions of the a— j3, fi—j, and a-7 in- 
terfaces, respectively. Finally, r is the line tension, i.e., 
the excess free energy per unit length of the three-phase- 
contact line we are interested in. Strictly speaking the 
line tension defined via Eqs. (|5.1| and (|5.2p is not iden- 
tical to its definition via Eq. (|1.1|) . because in Eqs. (|5.1j) 
and (|5.2p we first keep the subleading terms and treat 
them as subleading contributions to r, whereas in Eq. 
(jl.ip we isolate and drop these terms from the outset. 
However, if in the following we speak about the line ten- 
sion itself, as deduced from the decomposition of Q via 
Eqs. (|5.ip and ()5.2|) . we shall always drop the sublead- 
ing terms. In this sense the definitions in Eqs. (|5.ip and 
(|5.2p and that via Eq. (jl.ip do agree and we therefore do 
not introduce a different notation to distinguish between 
these two expressions of the line tension. 

For the lens the pressures in the a and 7 phases are 
equal, 

p a =p~ ( =p, (5.3) 

because the coexisting phases a and 7 are separated by a 
planar interface. For the drop the role of p 1 is played by 
— w 7 , i.e., the grand canonical free energy density of the 
solid while for the phase a we set p a = p. The pressure 
of the f3 phase deviates from p and is written as (see, c.f., 
the discussion following Eq. (|5.9[) ) 

p fJ = p + Ap. (5.4) 

The expression for the grand potential Q can be re- 
grouped as 

fi = + An , (5.5) 



11 



where 



for the lens, and 



Oo = —p V + A a a-, 



n = -p {v a + Vp) + + Aa an 



(5.6) 



(5.7) 



for the drop. Consequently one obtains 

AO = —l±.pVf3+A a p<j a fi+Af3 1 Of3 1 —'Kr 2 <7 a ~ ( J r2'n:rT (5.8) 
for the lens, and 

AO = -Ap Vp + A a pa a fi + (ap 1 - <T aj ) irr 2 + 2?rrT (5.9) 
for the drop. 

The equations given above hold for the two following sce- 
narios. In the first scenario, the lens or drop can exchange 
matter with the surrounding phase so that the chemical 
potentials in all phases are equal and the pressure in the 
(3 phase is determined by the chemical potentials. In the 
second scenario, nonvolatile lenses or drops are investi- 
gated and the volume of the lens or drop is prescribed. 
In this case the pressure in the (3 phase is not an inde- 
pendent thermodynamic variable but is determined by 
the amount of (3 phase. The term O (Eq. Ij5.6|) ) is in- 
dependent of the lens size and of the choice of dividing 
interfaces. This holds also for the drop (Eq. 15.7ft ) due to 
the transformation law given in Eq. (|2.3p which is valid 
if the a-"f and (3—j interfaces form a common plane, as is 
the case for the convention we have chosen. Since O as a 
physical quantity and Oo, as just argued, are independent 
of choices of dividing interfaces, also AO is independent 
of such choices. Thus in the following we can focus on 
AO. In order to extract the value of r from a known, 
e.g., calculated AO, or to infer information about pos- 
sible changes of r upon shifts of the dividing interfaces, 
at first it is necessary to specify all other quantities ap- 
pearing in Eqs. (|5.8|) and (|5.9|) . The geometrical quan- 
tities like volumes, areas, and lengths are defined once 
the dividing interfaces are chosen. The surface tension 
that of a planar interface between phases a and 
7. If both phases a and 7 are fluid it is known (see, 
e.g., Ref. 0|) that a aj is independent of the position of 
the dividing surface between these phases. If one of the 
phases, e.g., 7 is an inert solid phase, cr Q7 must be an- 
alyzed more closely; but we postpone this discussion to 
a later subsection (see, c.f., Subsec. 5.3.1) in which ses- 
sile drops are discussed separately. The quantities a a p 
and cr ) 3 7 are surface tensions of curved surfaces - again 
the case of a drop requires a separate discussion - and 
thus they depend on the corresponding radii of curvature. 



We indicate this dependence a a p(Ri), cr/3 7 (i?2) explicitly 
whenever it is necessary in order to avoid confusion. As 
already discussed above we postulate that the interfacial 
tensions appearing in Eqs. (|5.8|) and (|5.9|) have the same 
properties as their spherically closed counterparts, such 
as 



c Q /?(i?i) = a- af3 (oo) 1 



2^ 
Ri 



.It. 



(5.10) 



were <5^g is the Tolmann length for the a- (3 interface; 
subleading terms in 1/R\ are not considered here. An 
analogous equation holds for o- J g 7 (i? 2 ). Although in the 
following these expressions will not be be used explicitly 
it will be always understood that surface tensions depend 
on the radius of curvature in such a way (concerning the 
physical radius, see the corresponding remarks below). A 
related property of the pressure differen ce Ap is describ ed 
by the generalized Laplace equation 0, 1 1 19L J 



1201.1121 1 



Ap = 



2a af }(Ri) 

Ri 
2a^(R 2 ) 

R2 



dcr a 



di?i 



di? 2 



(5.11) 



where the terms in square brackets are quantities termed 
by Rowlinson and Widom Q notional derivatives. Such 
a term multiplied by a small [dRi] approximates the 
change in surface tension upon increasing the radius of 
the dividing surface by that value [dRi] without chang- 
ing the physical system, i.e., the density distributions 
and all thermodynamic variables remain fixed whereas 
the description of the system has been changed by shift- 
ing the dividing surface. In the following square brack- 
ets will be always used in order to characterize notional 
changes, e.g., [dRi] denotes the notional change of the 
radius whereas we would write dRi if we speak about 
a physical change of the radius at a fixed convention of 
choosing the dividing surface (interface). The quantity 
depends on the choice of the dividing surface. The 



dRi 



so-called surface of tension is defined as that dividing sur- 
face for which 



dRi 



is zero. For the equimolar dividing 
surface the notional derivative coincides with the deriva- 
tive of the surface tension with respect to the physical 
drop size, i.e., the derivative of Eq. (|5.10[) . The second 
term on the rhs of Eq. (|5.1ip renders the rhs invariant 
with respect to changes of the dividing surface. This 
must be the case because Ap defined as the pressure dif- 
ference between two bulk phases at given thermodynamic 
conditions is a measurable and therefore invariant quan- 
tity. Below, notional derivatives will be characterized by 
putting them into square brackets. 
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(2) 




(i) 


2a$„ [dRi] 


di?,_ 











A property that will also be used extensively in the 
following is the fact, that <j^ v (Ri) (i e {1,2}, £ 
{a,/?, 7}) is independent of the particular choice of the 
dividing surface up to and including the order 1/Ri- At 
some occasions we shall use the relation 



+ s.l.t. (5.12) 



between notional derivatives of surface tensions for two 
differently chosen dividing interfaces. The two dividing 
interfaces are denoted by the superscripts (1) and (2) 
and their radii are related via = [Ri]^ + [dRi]; 

[dRi] could be also understood as a differential. It could 
be also finite but then we assume that it is small com- 
pared with Ri. In Eq. (|5.12|l we write Ri and do not 
introduce or [i?;]'" 2 '' because using one or the other 

convention would lead to expressions which only differ in 
subleading terms which arc neglected anyway. Equation 
([5TT2"|) is valid if the values of both [R t ] {1) and [R t ] {2) are 
close to that of the radius corresponding to the surface 
of tension and it is a direct consequence of well-known 
res ults obtained for spherical interfaces (see, e.g., Refs. 
[7I ll2Clj |). The neglected s.l.t. contain terms of the order 
(1/Ri) 3 or higher and in addition they contain [dRi] or 
differences between and the radius corresponding 

to the surface of tension raised to the second or higher 
power. Most of the following conclusions, however, do 
not make explicit use of Eq. (|5.12[) . but only of the fact 



that 



~dR~ 



(2) 



(1) 



dcr£„ 




dcr ?l/ 


dR t 




dR t 



is of higher order than I/R4. 

We also use (in the second part of our reasoning) the 
following relation - known for spherical drops [J 1 120| 
- between the stiffness against changes of the radius of 
curvature at given thermodynamic conditions (i.e., fixed 
temperature and chemical potentials), denoted here as 
1 , and the notional derivative of the surface tension: 



(5.13) 



In the following we discuss lenses and drops. Although 
the arguments proceed similarly in both cases, we discuss 
them separately because we want to address certain spe- 
cial problems connected exclusively with the solid-fluid 
interfaces. 



B. The lens 

First, we investigate the consequences of the require- 
ment that the grand canonical potential of the total 
system is independent of particularly chosen dividing 



interfaces. Secondly, we derive - from a variational 
procedure at fixed lens volume - the equations yielding 
the contact angles as a function of the radius r char- 
acterizing the lens size. The variational problem is set 
up in a way which guarantees that the contact angles 
transform correctly upon notional shifts of the Gibbs 
dividing interfaces. The results following from these two 
procedures are then compared. 



1. Notional variation of the grand canonical potential 

Notional variations (indicated by square brackets) 
leave the grand canonical potential unchanged which in 
"differential" form is written as 



[dAfi] = 



(5.14) 



with AQ defined in Eq. (J5lj). The left-hand side of Eq. 
(|5.14|) gives the notional change of Ail upon notional 
shifts of the dividing interfaces away from some given - 
but essentially arbitrary - set of dividing interfaces, in 
the limit of very small notional shifts, characterized by, 
e.g., [dRi] (i = 1,2) (see below). In the above equation 
all notional variations of geometrical quantities are then 
expressed in terms of the notional variations [dRi] (i = 
1, 2) of the two radii of curvature Ri and R2 such as, for 
example, 



[dVfs] = 



dVp 



di?i 



[dRi 



dVp 
dR 2 



[dR 2 



(5.15) 



Similar equations hold for the notional variations of inter- 
facial areas [cLAq.^], [dAg-y], [d(7rr 2 )], and for the length 
of the three-phase-contact line [d(27rr)] as well as for 
the surface and line tensions. The notional derivatives 



der Q 



and 



d<rg 7 
dR 2 



of the surface tensions are - as al- 
taken to be identical to those defined for 

(Em - 1533) ) 



dRi 

ready stated 

completely spherical drops (see Eqs 



addition we demand that 



da a 



dR 2 



= and 



da 3 



dR t 



In 

= 



because a notional change of the radius R\ of the a— 
f3 interface should not lead to notional changes of the 
/3~7 interfacial tension, i.e., of the interface on the op- 
posite side of the lens, and vice versa (see Fig. 4). Wc 
further introduce two notional derivatives of the line ten- 



sion: 



dr 



diii 



and 



dR, 2 



These are actually defined by Eq. 
(|5.14p because all other quantities in this equation are ei- 
ther given by geometry or fixed via the set of definitions 
given above. Certain properties of these new quantities 
will follow from the analysis given below. 
The pressure difference Ap is related to the surface ten- 
sions and their notional derivatives via Eq. (|5.11|) . Since 
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Ap is unique the right hand sides of these two equations 
in Eq. (|5. 1 1[) must be equal. (The validity of the Laplace 
equation for even small lenses has been checked by sim- 
ulations, e.g., in Ref. |l5d |.) After using the geometrical 
relations R± = rj sin a and R 2 = r/sin.7 we obtain a 
first equation relating a, 7 and r: 



a a p sin a - 



dRi 



apy sin 7 + ■ 



di? 2 



Because the notional changes [di?i] and [6.R2] are inde- 
pendent, the prefactors of [di?i] and [di?2], obtained af- 
ter expressing Eq. (|5.14[) in terms of these variables, 
must both vanish. In this way two further equations 
are obtained leading to three equations in total. On 
the other hand only two equations are required in or- 
der to determine the contact angles for a given physical 
lens size, because the geometric parameters a, /?, and r 
are already related through a prescribed volume of the (3 
phase. Therefore the two aforementioned equations fol- 
lowing from Eq. (|5. 14|) have to be identical. This leads 
to the following consistency condition: 



cos 7 



" dr " 




" dr ' 




dcr / 3 7 






— cos a 




= i{ 




dR 2 


di?! 




di? 2 





da. 



a0 



di?i 



(5.17) 



The remaining third equation can be written in a sym- 
metrized form, 



a a p cos a + <7/3 7 cos 7 + a aj = 

r {sin 7 cos a — sin a cos 7} 
4 cos a cos 7 
{sin 7 cos a + sin a cos 7} 



da 



/3 7 



2 cos a cos 7 



cos 7 



di? 2 
dr " 



da. 



a/3 



di?i 



di? 2 



dr 



di?! 



(5.18) 



or, using the consistency condition given in Eq. (|5.17p . 
it can be put into the following form: 



a a fj cos a+a ) 3 7 cos j+a c 



-+sin7 



dr 



di? 2 



+sina 



dr 



dR 1 



(5.19) 

Since the contact angles a and 7 are natural variables 
characterizing the three-phase-contact region we express 
notional derivatives of r with respect to the radii in terms 
of notional derivatives with respect to the contact an- 
gles: 



dr 


= [ — 1 


da 


+ 


dr" 




d7 


, where 


da 


dR 4 


LdaJ 


dRi 




d7 




dRi 


dRi 



and 



d7 



dRi 



describe notional changes of a and 7, respec- 
tively, upon notional changes of Ri. For given centers M\ 
and M 2l the descriptions in terms of the pairs (a, 7) or 
{Rl,R>}) are equivalent. However, we do not transform 
the notional derivatives of the surface tensions, because 
the radii of curvature Ri are the natural variables of the 
curved interfaces. In the new variables Eqs. (I5.17|) and 



(5.16) (pTT9|) take the form 



"dr" 


. 2 


'dr" 


da 


sin a — 




d7 



sin 7 = 



and 

a a p cos a + a py cos 7 + a ai 
sin 7 cos 7 



da 



/37 



dR 2 



da. 



a 1 J 



dR 1 



t sin a cos a 
- H 

r r 



(5.20) 



dr 



dr 

d7 



da 
(5.21) 



The next question is whether the line tension r de- 
pends on the choice of the dividing interfaces. In order 
to answer it we compare Ail (Eq. (15.81) ) evaluated for 
two different dividing interfaces and find r to be indepen- 
dent of the choices of the dividing interfaces within the 
leading order. More precisely this statement says that 
notional shifts of the dividing interfaces change r only 
by contributions which decrease as 1/r with increasing 
r or faster. (We recall our remarks after Eqs. (|5.ip and 
(|5.2p saying that r as defined by Eqs. (|5.ip and (|5.2[) may 
contain subleading terms.) The leading contribution to r 
which is independent of r, is independent of the choices 
of the dividing interfaces. Although the subleading con- 
tributions to r itself may be neglected in the term r/r in 
Eq. (|5.19|) if we want to keep only terms up to the order 
1/r, it is not permissible to neglect the terms contain- 

If these derivatives 



dRi 



ing the notional derivatives 
are of the order 1/r - in fact we shall see below that 
notional shifts of the dividing interfaces lead to changes 
which are of that order - the pertinent terms 



LLL [dRi 

in Eq. (|5.19|) are of the order of 1/r, i.e., they are of 
the same order as r/r with only the leading term in r 
kept. In proving the aforementioned results concerning 
the behavior of r under notional shifts of the dividing 
interfaces, we have used the fact that surface tensions 
are independent of the chosen dividing interfaces up to 
order 1/R. Furthermore we have used the generalized 
Laplace equation (|5.1ip and also the fact that the differ- 

\d<7 1 (2) [da I (1) . 

ence -g^r — ~Kr l ^ s a correction of higher order 
than we are interested in. Moreover, the contact angles 
and other geometrical quantities for a choice (2) for the 
dividing interfaces have to be expressed in terms of the 
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corresponding quantities characterizing choice (1). This 
has been carried out up to the order giving rise to con- 
tributions of the order of the line-tension contribution. 
Next we investigate the transformation behavior of the 



quantities 



dr 



dRi 



due to notional shifts 



of the dividing interfaces. For this purpose we consider 
the equations for the contact angles (Eqs. (|5.16l - |5.21[) ) 
for two different sets (1) and (2) of dividing interfaces. 
The geometrical quantities, say a^ 2 \ 7^ 2 \ and r^ 2 \ in 
equations valid for the set (2) are then expressed in terms 



of a^ x \ 7^ 1 - ) , and up to the required order. The same 
is done for the surface tensions and their notional deriva- 
tives. We also use the above result that r is independent 
of the choice of dividing interfaces to leading order. From 
comparison of the two systems of equations (Eqs. (|5.16|) . 
(|5.17p . and (|5. 19|i ) and the fact that both systems must 
yield the same relations between r ) tyl \ and A 1 ^ for a 
fixed physical system, by using Eq. (|5.12|) we obtain the 
following transformation rules: 



" dr " 


(2) 


' dr " 


cLRi 




di?i 


" dr " 


(2) 


' dr " 


di? 2 




di? 2 



(1) 



a a p sin a cos(a + 7) 
r sin(a + 7) 

o~ a p sin a 
r sin(a + 7) 



dRi] 



r sm(a + 7) 

crg 7 sin 7 cos(a + 7) 
r sin(a + 7) 



[dR 2 ] , 
[dR 2 ] ■ 



(5.22) 



Alternatively, this may be translated into the notional 
derivatives of r with respect to the contact angles: 



"dr" 


(2) 


"dr" 


(i) 


da 




da 




"dr" 


(2) 


"dr" 


(i) 


d7 




d7 





= -CT Q(3 [di?i] 
= [dR 2 ] 



(5.23) 



Equations (|5.22|) and (|5.23|) show that the notional 
derivatives of r depend on where the dividing interfaces 
are located. This is completely analogous to what is 
known for the notional derivative of the surface tension 
for closed spherical interfaces. 

2. Variational treatment with constraint of fixed volume 

A common procedure aimed at obtaining equations for 
the contact angles consists of minimization of the grand 
canonical potential AQ with the constraint of fixed vol- 
ume Vp (compare Eq. ((| 



Afi = A a pa a p(Rx) + Ap 1 ijp 1 (R 2 ) - irr 2 a 
+27rrr(a, 7, r) + XVp 



(5.24) 



In Eq. (|5.24[) no bulk contributions appear because the 
volume Vp, the bulk pressures p a and pp, and thus the 
pressure difference Ap — pp — p a are kept fixed due to 
the fixed thermodynamic conditions specified by the tem- 
perature and the chemical potentials. The fixed volume 
condition is implemented via the last term in Eq. (|5.24p 



containing the Lagrange muliplier A. The volume Vp is 
the one enclosed by the dividing interfaces. The indepen- 
dent variables of the variational problem are a, 7, and r 
(f3 is determined via a + 7 + [3 = 2n) . In order to calcu- 
late the variation dAfi resulting from the variations of a, 
7, and r we have to introduce - in addition to the inter- 
facial and the line tensions - new material parameters. 
These are the stiffness constants -gjf^ | and \ of the 
interfaces describing the cost in interfacial free energy 
resulting from changes in the radii of curvature as well 
as the stiffnesses ^ \ , ^ | , and ^ | of the three-phase- 
contact line against changes of the contact angles a and 
7, and the radius of curvature r, respectively. To be ex- 
plicit: in the variation we introduce also contributions 
A a p-^^-\^rdr describing the cost in free energy associ- 
ated with a variation of the curvature radius i?i , caused 
by a change of r by dr, without changing the interfacial 
area, or 27rr^-|dr etc.. The symbol | indicates that the 
stiffnesses are given by the respective derivatives evalu- 
ated at fixed thermodynamic conditions. 
The stiffness constants of the interfaces, e.g., -^j-, — 
should not be confused with the derivatives of surface 
tensions with respect to the physical radius. For exam- 
Gaf> does not coincide with the derivative of a, 



d<T a /3 I 



d_Ri 



ap 



pie, 

given in Eq. (|5.10[) with respect to R\ because drops 
of different physical sizes in an unstable or constrained 
equilibrium with their environment do not correspond to 
the same thermodynamic conditions given by tempera- 
ture and chemical potentials. The stiffness constants of 
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the interfaces may be expressed via Eq. (|5. 13[) in terms 
of the notional derivatives. This relation shows that the 
stiffness constant of an interface vanishes if the so-called 
surface of tension is chosen as the Gibbs dividing surface. 
In the same way the stiffness constants attributed to the 
contact line should not be confused with similarly look- 
ing derivatives with respect to an implicit dependence of 
the line tension on contact angles. These implicit depen- 
dences reflect changes in the thermodynamic: conditions. 
The meaning of the stiffness constants is completely dif- 
ferent. In setting up the variational principle we explore 
constrained equilibrium configurations in the neighbour- 
hood of the equilibrium configuration in order to find the 
equilibrium by minimizing the free energy at fixed ther- 



modynamic conditions. Surface tensions, the line ten- 
sion, and the various stiffness constants describe the costs 
in free energy due to virtual displacements of interfaces 
away from their equilibrium shape. (In the summary op- 
erational procedures will be discussed how to determine 
stiffness constants theoretically.) 
The equation 



, . - dAfl , dAft , dAn 
dAQ = — — dr H — da H — 

dr da d7 



d 7 = (5.25) 



leads to three equations because a, 7, and r can be varied 
independently. They have the following form: 



2a, 



a/3 I 



da, 



a/3 



di?i 



(1 — cos a) f do - 



(1 + cos a) dRi 



?af3 



2 sin a dr 



da 



(5.26) 



and 



2fT/3 7 sin 7 dap-y (1 — COS7) J dayj 7 



a a /3 cos a + <7/3 7 cos 7 + c Q1 



dr 
dr 



di? 2 



cos a sin a dr 



2 sin 7 dr 



(1 + COS7) [ di? 2 

cos 7 sin 7 dr 



r da 
r(l — cos a) f da a p 



sin a 



di?i 



r d7 
2 sin 2 a dr 



sin 7 



do 



da 



r(l — cos 7) f der 



T /37 



di? 2 



(5.27) 



2 sin 2 7 dr 



d" 



(5.28) 



The ensuing equality of the right hand sides of Eqs. 
(|5.26|) and (|5.27|) leads to one of two equations relating 
the variables a, 7, and r. The second relation is pro- 
vided by Eq. (|5.28|) . Finally a third equation expresses 
the fixed volume Vp of the lens, 



1 3 J (1 + cos a) (2 — cos a) 
3 [ (1 — cos a) sin a 
(1 + cos7)(2 — COS7) 
(1 — cos 7) sin 7 



(5.29) 



and thus finally determines a, 7, r, and the Lagrange 
parameter A via Eqs. (|5.26|) or (|5.27|) . 

We now compare two such sets of equations, each for a 
different choice of the dividing interfaces in order to find 
relations between the stiffnesses of the line for two differ- 
ent choices for the dividing interfaces. Since the physical 
object is kept fixed, the relations between one set of vari- 



ables (a' 2 ), 7^, r^ 2 )) and another set {a^\ 7W, r^ 1 ') 
follow from the definition of the lens by two intersecting 
spheres (Mi,R{'' ) and (M 2 ,i? 2 ^) as described in Sec. 5 
and simple geometrical considerations. The volume Vp 
attributed to the lens depends on the choice of the divid- 
ing interfaces, i.e., two different volumes ^ Vp are 
assigned to the same physical object. The relation be- 
tween Vp and Vp is known once the notional changes 
of the radii R\ and R2, fixing the relative positions of the 
dividing interfaces, are given. Therefore Eq. (|5.29[) docs 
not contain any new information that could be used in 
order to relate the stiffnesses which have to be attributed 
to the line for two different choices for the dividing in- 
terfaces. Thus in the following, we have to consider only 
two pairs of equations: first, the equation resulting from 
the equality of the right hand sides of Eqs. (|5.26|) and 
(|5.27p . and secondly Eq. (|5.28[) . one pair for each choice 
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of the dividing interfaces. The parameters 

a (2) ;7 (2) 

, and 

r*- 2 ) can be expressed in terms of a^> , 7^, and and 
of the notional changes of the radii [dRi] — R^ — r\ , 
i = 1,2. After expressing the pair of equations for con- 



vention (2) in terms of a 



(i). 7 (i) 



and r-W up to the order 



of the line tension term, the resulting equations may be 
compared directly with the pair of equations obtained 
if convention (1) was chosen from the beginning. It is 
clear that the two pairs of equations must be equivalent 
since they describe the same physical system in terms of 
the same variables. In that comparison we further use 
the relation (|5.13|) between the stiffnesses of the inter- 
faces against changes of the radii of curvature and the 
related notional derivatives. In addition we use that the 
mentioned notional derivatives written for two different 
dividing interfaces are related via Eq. (|5.12| which in 
the chosen variables read (up to the relevant order it is 
not necessary to distinguish between and a*- 2 -* etc. 
and thus on the right hand side we omit the superscripts 
distinguishing between these conventions): 



s.l.t. 
(5.30) 



da a/ 3 


(2) 


da a p 


(1) 2cr Q(3 sin 2 a[di?i] 


dRi 




dRi 


r 2 



and 



da 7 


(2) 


derg 7 


(1) _ 2^ 7 sin 2 7 [dE 2 ] 


dR 2 




di? 2 _ 





+ s.l.t. 

(5.31) 

As a result of the comparison and the requirement of 
equivalence of the two pairs of equations discussed above 
(rhs of Eq. ([535)) equals rhs of Eq. ([537)1 . and Eq. 
(I5.28P ) we obtain two coupled equations for the follow- 
ing three quantities: (f^ 2 > - - 

and (gr|^ — 37 1^) • These equations have a manifold 
of solutions. If we pick a particular solution with 



dr | (2) 
dr ' 



dr ,(i) 
dr ' 



= 



(5.32) 



the two remaining quantities are given by 
= -a af} [dRi] 



dr 1(2) 
da I 



dr I 
da I 



- 2 I 


da af 3 
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da a 


2 sin 2 a I 






dfli 



±L|(2) _ dr 1(1) 
d7 I d7 I 



= -CT / 3 7 [dR 2 ] 



- 2 ( 


dep., 


(2) 


dap-, 


2 sin 2 7 I 


dR 2 




di? 2 



(1) 



(1) 



(5.33) 

For the special choice in Eq. ([5.32]) leading to Eq. (|5.33j) 
the Lagrange multiplier A (see Eqs. (|5.26[) and (|5.27|) ) 
becomes independent of the chosen dividing interfaces 
and therefore it allows for a physical interpretation. The 
conditions (|5.33|) are obviously fulfilled if the following 



relations between stiffness constants of the line and those 
of the interfaces hold: 



dr 
do 



da a /3 



2 sin 2 a dRi 



2 sin 2 a 



da, 



a/3 



dRi 



and 



dr 

do 



da 



0j 



2 sin 7 dR 2 



da 



dR 2 



(5.34) 



(5.35) 



2 sin^ 7 

Additional terms which are independent of the choice of 
dividing interfaces could be added on the right hand sides 
of Eqs. (|5.34[) and (|5.35[) without violating the conditions 
(|5.33[) . However, keeping such terms would not lead to 
more general results, but to equations which are more 
complicated than those presented below. As we shall see, 
the special choices of the right hand sides of Eqs. (|5.34|) 
and (|5.35[) lead to equilibrium conditions for the contact 
angles which agree with Eqs. (|5.16[) and (|5.21|) which 
have been derived above from the principle of invariancc 
of fl under notional changes. This agreement is achieved 
for the particular choice af| (2) - = 0, but is 

still left undetermined by consistency requirements alone. 
Insertion of Eqs. ([534]) and ([535)1 into Eqs. ([536]) and 
([537)1 leads to 



A 



and 



•la, 



a/3 I 



da, 



a/3 



di?i 



2a, 



a/3 



Rl 



da, 



a/3 



dRi 



X 



2a 



/3 7 



I? 2 



dap j 



dR 2 



(5.36) 



(5.37) 



In other words, A equals minus the Laplace pressure ex- 
pressed via the generalized Laplace equation (compare 
Eq. (|5 . 1 If) ^ which supports the choice made in Eqs. 
(|5.34p and (|5.35[) above. Accordingly, the equations re- 
lating the contact angles to the radius r read (recall that 
this holds only for choices which are in accordance with 
Eq. ([532)1 as well as Eqs. ([5T34| and ([535)) V. 



a a f3 sin a 
and 



da, 



a /3 



di?i 



CT3 7 sin 7 + 



da 



Pi 



dR 2 



(5.38) 



r dr 

a a a cos a + arjy cos 7 + a ay = — h — 

r dr 



sin a cos a dr 



da 



sin 7 cos 7 dr 



(5.39) 



Comparison of Eq. ([539)1 with Eq. ([531)1 (Eqs. ([535)1 
and (|5.16p are identical anyway) renders the relation 



dr 
dr 



sin a cos a dr 



da 



sin 7 cos 7 dr 



sin a cos a 



dr 



da 



r d7 
sin 7 cos 7 
r 



dr 

d7 



(5.40) 
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Adopting this relation for the surface of tension (s) and 
using the relations in Eqs. (|5.34[) , (|5.35[) , and (|5 . 32[) and 

00 



noting further that -jfrf- = = and thus 

= o = dz|( s ) we find (note that because of Eq. 
(|5.32[) is independent of the choice of the dividing 
interfaces) 



dr 
dr 



dr 
dr 



sin a cos a 



dr 



lW 



da 



sin 7 cos 7 



dr 
d7 
(5.41 



00 



On the rhs of Eq. (|5.4ip superscripts w are omitted for 
a, 7, and r because differences between the values of a, 
7, and r for different dividing interfaces give rise only to 
higher order corrections. Reinserting Eq. (|5.41[) into Eq. 
(I5.40p we further find 



dr 




"dr" 




'dr" 


da 




da 




da 



00 



and 



dr 




"dr" 




'dr" 


d7 




d7 




d7 



(5.42) 



(5.43) 



Of course, due to these relations the transformation be- 
havior for the stiffness constants of the line tension (Eqs. 
(|5.32[) and (|5.33p ) is consistent with that for the notional 
derivatives (Eq. (|5.23p ). It should also be noted that 
= and = . 

da 1 C17 



C. The drop 

1. Notional variation of the grand canonical potential 

In the case of a drop our analysis also starts from 
the requirement that the grand potential must be invari- 
ant with respect to notional changes. They consist of a 
change of the radius of the a-(3 interface by [dR] and 
a common shift of the oe-j and /3~7 interfaces by [dh]. 
From the requirement [dAf2] = we obtain two equations 
( [dR] and [dh] can be chosen independently) relating the 
contact angle 9 to the radius r. However, these equations 
have to be identical because already a single equation is 
sufficient to determine the relation between 9 and r. This 
leads to the consistency relation 



"dr" 




"dr" 
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da a p 




+ cos 9 






dh 




dR 


= ~2 


dR 



and the following equation for the contact angle: 

dr' 



a a /3 cos 9 + (a Pi ~ Va*i) = sin ( 



dR 



(5.44) 



(5.45) 



Similarly as for the lens we have introduced notional 
derivatives [gr] an d [gjj] • We also make use of a no- 
of the spherically shaped a-(3 in- 
we do not take into ac- 



tional derivative 
terface but 



■ at the same time 
count a notional derivative 



d(g y 3-y- 



dh 



of the difference 



in the interfacial tensions of the two planar substrate- 
liquid (7-/3) and substrate-vapor (7-a) interfaces. This 
calls for a comment. Here the quantity cr^ 7 — cr Q7 denotes 
the difference of two interfacial tensions corresponding to 
a situation in which the same substrate phase 7 remains 
in contact with either the a or the (3 phase. Both phases 
a and (3 are at the same pressure p. This quantity to- 
gether with the a— (3 surface tension at infinite radius de- 
fines — via the Young equation — the contact angle 9o 
of a macroscopicly large liquid drop (phase f3) in contact 
with its vapor (phase a) on top of the substrate 7. The 
angle 9q is an observable quantity and does not depend 
on the choice of the dividing surface. Since a a p does not 
depend on the choice of the dividing interfaces either, it 
follows that apy — a aj cannot depend on the chosen po- 
sition h of the common a-7 and (3-^f dividing interfaces. 
One arrives at the same conclusion by calculating the 
changes of cr^ 7 and a ai with respect to changes dh in 
the interface positions. In the difference cr^ 7 — a a ^ terms 
depending on dh drop out provided the pressure in the 
phases a and (3 is the same, and a common height 'above' 
the substrate is chosen for the a-7 and for the f3-j di- 
viding interfaces (see, e.g., the discussion at the end of 
Sect. 2). (Both these assumptions are implicitly con- 
tained in the Young equation.) Actually, further compli- 
cations would arise if, in an attempt to mimic the sit- 
uation in a drop more closely, we would carry out the 
calculations at different pressures, e.g., for pp = p a + Ap 
with the Laplace pressure Ap. A quantity (a^ — o ai )' 
defined in that way would depend on the choice of the 
dividing interfaces and additional terms of the same or- 
der as those introduced by the line tension would ap- 
pear in the final equations. Furthermore the difference 
(cr^ 7 — er Q7 )' would not be related to 9q via the standard 
Young equation and thus could not be determined via 
an independent experiment of measuring the shape of a 
macroscopic drop. Therefore it does not appear to be 
useful to introduce quantities like (crp-y — c Q7 )' (see, how- 
ever, the discussion in Subsect. VI C). 
Equations (|5.44[) and (|5.45|) may be rewritten in terms of 
9 and r. In these variables the consistency equation and 
the modified Young equation take the following form: 



sm 



"dr" 


r 2 


da a p 


dff 
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(5.46) 
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and 



7 a p COS 9 + (cr / 3 7 - <7 a7 ) = 



dr 
dr 



dr 
d0 
(5.47) 



Since the 7— a and 7-/3 interfaces are planar, a notional 
derivative of r with respect to the angle 7 as introduced 
for the lens cannot be defined in the present case; instead 
is used. As in the case of the lens the relations be- 
tween ([di?], [dh]) and ([d9], [dr]) are unique and the 
same is true for the corresponding notional derivatives. 
Similarly as for the lens we find that r is independent of 
the choice of dividing interfaces. This statement is true in 
the same sense as the corresponding one discussed for the 
lens below Eq. (|5.21|) . As in the case of the lens, higher 
order corrections to r decreasing with increasing r like, 
e.g., 1/r — such corrections are possible if r is defined via 
Eq. (|5.2|) instead of Eq. (jl.ip - depend on the choice 
of the dividing interfaces. Correspondingly, the notional 
derivatives of r appearing in Eqs. (|5.46p and (|5.47[) also 
depend on the choice of dividing interfaces and the cor- 
responding terms are furthermore of the same order as 
the term r/r with only the leading term in r being kept. 
We now analyse the transformation of the notional 
derivatives upon notional shifts of the dividing interfaces. 
The same kind of arguments as for the lens leads to the 
following relations: 
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(5.48) 



Changing the set of independent variables from (h, R) to 
(r, 9) leads to the following form of the above transfor- 
mation laws: 
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u a j3 sin 9 
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[dh] . (5.49) 



2. Variational treatment with constraint of fixed volume 

Here the objective is to minimize the grand canonical 
potential Afl with the constraint of fixed volume V@: 

Af2 = A a pa a p(R) - nr (a aj - crg 7 ) + 2nrT(9, r) + XVp 

(5.50) 



in which no bulk terms are included since the Laplace 
pressure and the volume of the (3 phase are considered to 
be fixed; the Lagrange multiplier is denoted by A. We use 
9 and r as the variables describing the system, and we in- 



troduce the stiffness constants 



da a 



and 



d<r Q 



which are 



dr I d8 

then expressed in terms of the 'natural' stiffness constant 

dR 



with known properties. In addition, similarly as for 
the lens, we introduce the stiffness constants ^| and ^| 
describing the cost in free energy attributed to the three- 
phase-contact line resulting from variational changes of 
r and 9 at fixed thermodynamic conditions. In the fol- 
lowing we shall determine how these stiffness constants 
depend on the choice of dividing interfaces and how they 
are related to notional derivatives. 

If we impose the relation (for a justification see the dis- 
cussion following, c.f., Eq. (|5.54j) ) 



sm 



d9 



da a 



2 di? 



(5.51) 



as in the case of the lens, the Lagrange multiplier A ac- 
quires the meaning of the negative Laplace pressure: 



A 



-Ap 



2a, 



a/3 ! 



1 9 da, 



a/3 



di? 



2a, 



a0 



R 



da, 



a/3 



di? 

(5.52) 



Equation (|5.5ip corresponds to the relations in Eqs. 
(|5.34[) and (|5.35[) which were introduced in the case of 
the lens for similar reasons. We also remark that Eq. 
(|5.51[) is formally identical to the consistency relation in 
Eq. (|5.46|) . The second equation, obtained from mini- 
mizing Afl and due to the independence of the variables 
9 and r, has the following form (compare Eq. (|5.28p ): 

dr 
r dr 
(2 + cos 9) sin6>dr 



r sm 



a a p cos 9 -+ 
da a p 



(°>y 



^7 J 



(l-cos0) di? 



d9 



(5.53) 



(An equation which expresses the fixed drop volume V/3 = 
(l/3)7rr 3 [(l-cos6»)(2+cos6»)]/[(l+cos0) sin 9] in terms of 
r and 9 and which together with the other equations fixes 
r, 9, and the Lagrange parameter A, does not contain 
information that could be used in the following. The 
same has been found for the lens.) 

Similarly as for the lens, we deduce from a comparison of 
the two equations, which follow from Eq. (|5.53p for two 
differently chosen dividing interfaces, the transformation 
laws between the stiffness constants of the three-phase- 
contact line: 
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(5.54) 
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Actually, the transformation behavior of ^ 
fixed independently via Eq. (|5.51|) . 



da a 



dR 



is already 
dCT ^" and 



dR 



Eq. (|5.12p . There is no contradiction between these two 
independent requirements because of the formal equiva- 
lence of the transformation laws in Eqs. (|5.49|) for the 
notional derivatives and those for the stiffness constants 
in Eqs. (|5.54[) . and because Eqs. (|5.46p and (|5.5ip are 
formally identical. In other words, one is indeed free to 
fix ^ by Eq. (|5.51[) and thus to bestow the meaning 
of the negative Laplace pressure —Ap on the Lagrange 
multiplier A. 

Equation (I5.53|) can be simplified by using Eq. (|5.51|) : 



<7 a p COS 9 + (cr ( g 7 — (T Q7 ) 



dr 
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sin f cos ( 



d9 

(5.55) 



From a comparison of Eqs. (|5.55|) and (|5.47[) we obtain 

(5.56) 
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(Note that applying Eq. (|5.51|) to Eq. (|5.56|) in combi- 
nation with Eq. ([533]) yields ^| (s) = 0.) Again, these 
relations are compatible with the transformation laws for 
both the notional derivatives and the stiffness constants. 



VI. POSSIBILITY OF CHOOSING DIFFERENT 
DEFINITIONS FOR LINE TENSIONS 



At this point it is due to comment on the seeming con- 
tradiction between our statement that r is independent 
of the choice of dividing interfaces and Eq. (|3.3|) , follow- 
ing from a discussion of a straight three-phase-contact 
line at a wedge shaped fluid volume, which states the op- 
posite and tells how r should change upon shifting the 
fluid-solid dividing interface. A question related to that 
issue is whether one can use the line tensions of a gas- 
liquid-solid contact calculated from a certain microscopic 
theory in combination with a decomposition scheme in a 
model system like that considered in Sect. Ill (i.e., for a 



fluid wedge geometry, see, e.g., Refs. [19J, l20[), and use 
them in the formulae derived above for the drop. 
In order to answer these questions we reanalyze again the 
particular wedge geometry investigated already in Sect. 
Ill in order to find out whether implicit assumptions in 
the analysis given there are at variance with prescriptions 
used in Sect. V for the drop. We also resume our discus- 
sion of the drop system in order to see whether sensible 
alternative definitions of the line tension in that system 
do exist and whether the contradiction mentioned above 



can be resolved by using an alternative definition of r. 
If such an alternative exists the question about the re- 
lation between the two different line tensions has to be 
answered. 



A. Decomposition and reassembly of a system: 
difficulties with edge terms 

In Sect. Ill we have calculated the changes of the areas 
of the liquid-gas and of the solid-fluid interfaces due to 
a notional shift of the position of the solid-fluid interface 
within the box with rectangular cross section, as shown 
in Fig. 2, implicitly assuming that these area changes 
are representative for an entire system from which Fig. 
2 actually shows only a small part. In this procedure a 
possible source of errors might be hidden. We therefore 
discuss now the underlying general procedure implicitly 
used in Sect. III. First, the total system is decomposed 
into two subsystems I and II such that subsystem I con- 
tains the three-phase-contact line as the object of inter- 
est. On the other hand subsystem I is surrounded by a 
subsystem II which itself does not contain any actual line 
contribution. (In the case of straight contact lines this 
requires to use periodic boundary conditions in one direc- 
tion; this is of no relevance for the following discussion.) 
The contribution of subsystem I to the grand potential 
is then analyzed in terms of volume, interfacial, and line 
contributions ignoring unphysical edge terms. The same 
is done for subsystem II except that no actual line con- 
tribution appears in the decomposition. At the end all 
contributions from both subsystems are put together in 
order to calculate the grand potential of the total system. 
It is expected that the grand potential (and the decom- 
position into volume, interfacial, and line contributions) 
calculated in that way is equal to what one obtains if one 
decomposes the grand potential of the whole system into 
volume, interfacial, and line contributions and then adds 
up these contributions without first partitioning the sys- 
tem. In what follows we discuss why it may be sometimes 
misleading to ignore the unphysical edge terms. 

The decomposition introduced above is not unique. 
Two examples of such decompositions are indicated in 
Fig. [6] by two boxes of different shapes enclosing and 
defining the respective subsystem I. The structure is as- 
sumed to be translationally invariant in the direction per- 
pendicular to the plane shown in Fig. [6l The cross sec- 
tion of the box has the shape of a rectangle for one of the 
chosen decompositions and in the other case it has the 
shape of a parallelogram such that two of the faces of the 
box are perpendicular to the liquid-gas interface. The 
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embedding system II is not specified in any detail except 
that it is understood that the structures found within the 
box, e.g., the liquid-gas interface, extend beyond the box 
boundaries into the embedding system II. For that rea- 
son interfacial energies at the boundaries of the boxes do 
not arise. (This statement also comprises the solid-liquid 
and the solid-gas interfaces since the respective bound- 
ary of the box can be placed anywhere inside the solid.) 
There are also no actual line contributions to the free en- 
ergy from the edges of the box or from the line at which 
the liquid-gas interface intersects a face of the box. The 
environments of these lines are not different from bulk 
or from a corresponding region in an infinitely extended 
interface. Inspite of the absence of actual line contri- 
butions, artificial non-physical line contributions at the 
lines just mentioned have to be included in the expres- 
sions for the grand canonical potentials characterizing 
the subsystems. That there is a need to introduce these 
artificial line terms becomes very plausible if one looks 
at the two possible decompositions of a system into two 
subsystems shown in Fig. [6] Evidently interfaces are 
cut in different ways by the faces defining subsystem I 
for the two considered options. Therefore, contributions 
to the grand canonical potential stemming from certain 
parts of the interfaces are either attributed to subsystem 
I or subsystem II depending on how the decomposition 
is done. Only if one includes the artificial line terms in 
the analysis one can keep track of such traces of the de- 
composition; otherwise one runs into inconsistencies. 

Before we proceed with this analysis, which ultimately 
aims at finding out whether and how the true line tension 
may be determined from microscopic calculations for the 
most simple geometries, we continue with the additive 
decomposition of the grand canonical potential f2 of a 
system 

n = n 1 + Q u , (6.1) 

where f2i and Qu are the contributions to the grand 
canonical potential attributed to subsystem I and II, 
respectively, which, e.g., in the spirit of density func- 
tional theory can be thought of to be calculated from 
the density distributions characterizing the total sys- 
tem. Already this decomposition is problematic and not 
unique although obvious choices do exist. For instance, 
if we base our theoretical description on a type of model 
in which the grand canonical potential is expressed in 
terms of a local functional of the number density (of, 
e.g., the Ginzburg-Landau type), f^i is obtained by in- 
tegrating the local density of the grand canonical poten- 
tial over the volume of the box defining the subsystem 
I. The local density of the grand canonical potential is 



obtained from particle-density distributions which may 
be obtained from calculations restricted to the interior 
of the box. The boundary conditions have to correspond 
to a seamless continuation of the structures in the inte- 
rior of subsystem I into the surrounding subsystem II. 
The decomposition procedure becomes less obvious if a 
non-local density-functional description is used. In that 
case interactions across the boundaries of the subsystems 
occur. In order to calculate the contribution to f2i com- 
ing from two-center integrals one might, e.g., restrict one 
of the integrations to the interior of the box, whereas 
the second integration runs over the total volume of the 
whole system. The value obtained for Q\ is independent 
of size and shape of the surrounding subsystem II pro- 
vided it is sufficiently big and the interactions between 
particles decay sufficiently rapidly as a function of their 
distance. Artificial interfacial contributions from the sur- 
faces of the box are eliminated automatically by chosing 
the described procedure. By carrying out the integra- 
tions in the way outlined above we basically mimic a 
local density of f2 allowing, however, for the embedding 
of the subsystem into a global system. 

In the sense of the previous paragraph f2 can be decom- 
posed into a sum of two contributions originating from 
two subsystems I and II. Since for this decomposition no 
use has been made of the concept of dividing interfaces 
between phases, in particular between solid and liquid or 
solid and gas, it is also clear that the two contributions 
r2j and On must be independent of choices for the divid- 
ing interfaces between solid and fluid. 
In the next step f2i and Qu are further decomposed into 
volume, interface, and line contributions, i.e., 

(h = Viuj + A\ g ai g + A^ct s i + ^ g cr sg + Lt + ^2 LiT^' 1 , 

(6.2) 

and 

flu = V u lo + Aftoig + Al\a sl + AHa se + £ L iT f' 11 , 

(6.3) 

where A\ and A\ g are the areas of the 1-g interfaces 
within subsystem I and II, respectively, A^ and A^ are 
those of the s-1 interfaces, and A\ g and A l ^ g are the areas 
of the s-g interfaces (1: liquid, g: gas, s: solid or wall). 
t is the line tension of the actual three-phase-contact 
line and thus it appears only in subsystem I. The terms 
t"*' 1 or Ti art ' n characterize the artificial line contribu- 
tions resulting from the decomposition into subsystems. 
The quantities Li are the lengths of the artificial lines, 
which are equal in both subsystems. Since Q = Oi + Qn 
and because the only actual line inhomogeneity, which is 
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FIG. 6: Two possibilities for cutting a subsystem I out of 
a total system containing a 1-g interface and a three-phase- 
contact line with the solid s. These two possibilities are in- 
dicated by two boxes (full and dotted lines) enclosing the 
1-g-s contact line. The structure is assumed to be transla- 
tionally invariant perpendicular to the plane of the figure. 
Two possible choices for the solid-fluid dividing interface are 
indicated by solid (1) and dashed (2) lines. The two dividing 
interfaces are separated by a distance Sh. The areas of the 
solid-liquid (solid-gas) interfaces within subsystem I change 
when the solid-fluid dividing interface is shifted. The values 
of these area changes depend on the shape of the box defining 
subsystem I. This dependence is indicated by the two hori- 
zontal double arrows. The short double arrow indicates the 
changes, due to the shift Sh, of the areas of the solid-gas 
and the solid-liquid interfaces within the box of rectangular 
cross section, the long double arrow the changes of these areas 
within the box with parallelogram-shaped cross section. 

present in the total system, is accounted for by the line 
tension r, one has 



art, I 
T- = — 



E T art, II 



(6.4) 



i.e., the artificial line contributions of the two subsystems 
cancel each other. 

Since f^i is independent of the choice of the dividing 
interfaces between phases, in particular of the solid-fluid 
dividing interface, the same arguments as used in Sec. 
3 can now be applied to the subsystem I and to f^i in 
order to study the influence of a shift by an amount Sh 
of the s-g and s-1 dividing interfaces. (We restrict the 
discussion to a structure that is translationally invariant 
in the direction parallel to the line, i.e., Li — L.) If the 
rectangular box indicated in Fig. [5] is chosen one obtains 



-1 (2) 



E T art, I 



1 



_art,I 



= <7i g sin OSh 
(6.5) 

for the difference of the expression in square brackets 
for the two choices (2) and (1) of the dividing interfaces 
(compare Fig. 2). The only but essential difference to 
Eq. (|3.3[) is that Eq. (|6.5[) gives the difference not only 
for t but for the sum of the actual and the artificial line 



energies. There seems to be a chance that the seem- 
ing contradiction between the result that r of the actual 
three-phase-contact line is independent of the choice of 
dividing interfaces and Eq. (|3.3p can be resolved by in- 
cluding the artificial line contributions and replacing Eq. 
([OJ) by Eq. ([53]) . Of course Eq. JB3J) is only valid 
for the decomposition employing the rectangular box. If 
instead of a rectangular box a box with a cross section 
of the shape of a parallelogram (as shown in Fig. [5]) is 
chosen Eq. (|6.5p is replaced by 



5> 
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(6.6) 

because the changes of the solid-liquid (solid-gas) inter- 
faces within subsystem I, resulting from the shift of the 
solid-fluid dividing interface, are computed in the way in- 
dicated by the long dotted double arrow in Fig. [5] and not 
in the way indicated by the short solid double arrow in 
Fig. [5] which is the correct way to compute those changes 
within the rectangular box. From the comparison of Eq. 
(|6.5|) and Eq. (|6.6|) it becomes obvious that in order 
to deduce r from a calulated expression for Or one has 
to subtract not only volume and interface contributions 
but also the artificial line contributions. Such artificial 
line contributions are generated if a box boundary cuts 
through the inhomogeneous structure of an interface in a 
non-adapted way. Within an interfacial region, which has 
a macroscopic extension in the lateral directions, the den- 
sity profiles of the constituents of the fluid only depend 
on the coordinate perpendicular to the interface. If the 
box boundary defining subsystem I cuts perpendicularly 
through the interface, the cut is adapted to the density 
profiles and no artificial line contribution is generated. 
If by contrast, a box boundary cuts through the inter- 
face with an angle deviating from 90°, a spatial region 
of columnar shape, filled with an inhomogeneous fluid, is 
generated above or below the interface, which no longer 
belongs to subsystem I although according to the interfa- 
cial area attributed to subsystem I it should contribute to 
the interfacial energy within that subsystem. Vice versa 
a volume of columnar shape is attributed to subsystem I 
although the fluid within that volume should contribute 
to the interfacial energy attributed to subsystem II. Since 
in general the free energy contributions from these two 
columnar spatial regions do not compensate one is left 
with an artificial line contribution. 

For that reason we discuss still another possibility (Fig. 
[7]) for cutting a subsystem I out of a total system con- 
taining a 1-g interface and a three-phase-contact line with 
the solid s. For this choice the box boundaries cut per- 
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pendicularly, i.e., in an adapted manner through both 
interfaces. This choice avoids the appearance of artifi- 
cial line terms at the cuts of the box boundaries through 
the liquid-gas or the solid-fluid interfaces. An additional 
line generated by this particular choice of the box bound- 
aries is placed such, that the system is completely homo- 
geneous in a big volume around that line. Therefore, no 
net artificial line contribution is associated with that line, 
too. Carrying out the same analysis as for the two other 
box shapes we again obtain Eq. (|6.5|) but the artificial 
line energies should now vanish leading back to Eq. (|3 . 3|) . 

After our previous analyses for different box shapes 
defining subsystem I one might come to the following sug- 
gestion to cope with artificial line contributions. For the 
rectangular box (which leads to Eq. (|6.5j) ) the only arti- 
ficial line contribution stems from the intersection of the 
liquid-gas interface with the box boundary. No artificial 
line contributions arise at the cuts of the box boundaries 
with the solid-fluid interface, because the box bound- 
aries cut perpendicularly through that interface. There 
are also no artificial line contributions from the edges of 
the box located in homogeneous regions of the fluid. If 
we now argue that the artificial line contribution does 
not change upon a notional shift of the solid-fluid inter- 
face because the corresponding line is far away from that 
interface, the terms related to the artificial line contri- 
bution would cancel in Eq. (|6.5j) and we would again 
be led back to Eq. (|3.3p . In the case of the parallelo- 
gram shaped box and Eq. (|6.6|) we might argue using the 
same type of arguments, that only the two line contribu- 
tions due to the intersection of the box boundaries with 
the solid-fluid interface do arise. But now one would 
be inclined to admit that these artificial line contribu- 
tions might change upon a a notional shift of the solid- 
fluid interface because the lines are spatially associated 
with the shifted interface. Therefore, one would not draw 
conclusions about notional changes of r from Eq. (|6.6p . 
However, one should not rely too heavily on such type 
of arguments. After all, notional shifts of dividing inter- 
faces lead just to changes in interfacial areas and thus 
to changes in interfacial contributions to the free energy 
which are then subtracted from the total free energy in 
order to define the line tension(s). From the viewpoint 
of a macroscopic theory it is not clear how these changes 
of interfacial free energies should be split up among the 
real and the artificial line contributions. 
In the present subsection we have discussed how a line 
tension r and its dependence on the choice of dividing 
interfaces can be determined unambigously from calcu- 
lations for a simple model system, i.e., a liquid wedge 
in contact with a gas phase on a solid substrate. In or- 



der to be able to proceed, one has to carry out the cal- 
culations within a finite subsystem (box) cut out from 
the unbounded model system. The seamless continua- 
tion into the embedding system has to be incorporated 
into the calculation by choosing proper boundary condi- 
tions and by taking into account certain interactions of 
the finite subsystem with the embedding system. In ad- 
dition, the box boundaries have to be choosen such, that 
they cut through interfaces perpendicularly and artificial 
box edges have to be placed into homogeneous regions of 
the total system (a possible choice of the box is shown 
in Fig. [7]). With all these precautions we find that the 
dependence of r on the choice of the solid-fluid dividing 
interface is definitely given by Eq. 1|3.3|) , 




FIG. 7: Another possibility for cutting a subsystem I out of 
a total system containing a 1-g interface and a three-phase- 
contact line with the solid s. The structure is assumed to be 
translationally invariant perpendicular to the plane of the fig- 
ure. Two possible choices for the solid-fluid dividing interface 
are indicated by solid (1) and dashed (2) lines. The two divid- 
ing interfaces are separated by a distance Sh. The areas of the 
solid-liquid (solid-gas) interfaces within subsystem I change 
(horizontal double arrow) if the solid-fluid dividing interface 
is shifted. These changes in areas are identical to those ob- 
tained if the box of rectangular cross section is chosen in Fig. 

m 



B. Finite containers filled with fluids 

After the discussions in the previous subsection we are 
still left with the contradiction between Eq. (|3.3[) . which 
tells us how r should depend on the choice of dividing 
interfaces, and our statement in Sect. V which says that 
r should be independent of that choice. We would like 
to stress that this latter statement is based on a well de- 
fined decomposition scheme of the grand potential of a 
drop on a substrate (or a liquid lens at a fluid interface) . 
In order to gain further insight into this problem we now 
discuss a finite fluid system enclosed in a container with 
solid walls. This way we avoid from the outset the neces- 
sity to cut out a finite subsystem and we stay away from 
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the danger to pick up unphysical edge or line contribu- 
tions. Our treatment of closed finite containers gives us 
also the opportunity to discuss a further aspect of solid- 
fluid systems, i.e., the curvature of the solid-fluid inter- 
face. We further chose the boundary conditions such that 
the liquid-vapor interface is planar and thus no compli- 
cations occur related to the Laplace pressure. In addition 
there is no curvature correction to the liquid-vapor sur- 
face tension. On the other hand if a container with only 
planar walls is partially filled with a liquid in contact 
with its vapor phase, in addition to the liquid-vapor- 
solid contact line further edge (line) contributions do ap- 
pear which scale with the same linear dimension of the 
container as the three-phase-contact line we are inter- 
ested in. Therefore, the different line (edge) contribu- 
tions cannot be separated unless these edge contributions 
have been determined independently from investigations 
of reference configurations not containing a liquid-vapor- 
solid contact line. As a reference configuration the same 
container but completely filled with either liquid or gas 
(vapor) is chosen. But then one encounters the problem 
that there is no unique description how to extract all ad- 
ditional edge contributions individually and in particular 
it is impossible to find the transformations upon notional 
shifts of dividing interfaces for each of these edge contri- 
butions individually. 

In order to facilitate to vary the length of the liquid- 
vapor-solid contact line indepently from the lengths of 
container edges we consider a biconical container com- 
posed of two identical but oppositely oriented right cones 
with circular base which are glued together base to base 
along the circumference of the base as shown in Fig. [8] 
All container walls are composed of the same solid ma- 
terial s. The cones are filled with liquid 1 up to a certain 
height, which is determined by the amount of liquid pro- 
vided, and the liquid is in contact with its vapor phase g. 
The opening angle of the cone is chosen to be 180° — 26 
such that the liquid-vapor interface is planar. The sys- 
tem contains two structural units characterized by lines. 
The first line is at the joint of the two cones and is termed 
Lcont (Fig- [£])• It has the shape of a circle with a radius 
denoted as r. In the following the related line (edge) 
tension will be called r con t. The second line, the one we 
are interested in, is the three-phase-contact line between 
solid, vapor, and liquid. It also has a circular shape with 
a radius denoted as r\ g . The related line tension will 
be called r. Again a parallel displacement of the divid- 
ing interfaces between solid and fluid by an amount 5 
is considered. This parallel displacement is indicated in 
Fig. [8] This set-up allows one to vary the length of 
the three-phase-contact line corresponding to r indepen- 




FIG. 8: Cut through a biconical container (circular base) with 
an opening angle 180° — 29. The apexes of the container are 
filled with the liquid phase (1) of a fluid up to the planar 
liquid-gas interfaces. The central part of the container is filled 
with the coexisting gas phase. The filling height is determined 
by a prescribed amount of liquid. Two possible choices for the 
dividing interface between solid (wall) and fluid are indicated 
by solid and dashed lines. The two possible dividing inter- 
faces are separated by a distance S. The two corresponding 
positions of the liquid-gas-solid contact lines denoted as L 
are indicated by small and large dots. At the joint between 
the two cones there is another circular line L cont formed. 



dently from the container dimensions. It is tempting to 
conclude that thus the line tension r and its dependence 
on the choice of dividing interfaces can be separated un- 
ambigously from the line contributions related to the ge- 
ometry of the container. However, since the container 
walls are curved, we also have to take into account cur- 
vature corrections to the solid-fluid interface tensions. 
As we shall see, these curvature corrections bring about 
contributions to the grand potential which scale with the 
length of the liquid-vapor-solid contact line. Therefore, 
for this geometry one also has to study at first reference 
configurations not containing a liquid-vapor-solid con- 
tact line, i.e., the bicone being homogeneously filled with 
either liquid or vapor, in order to determine how the cur- 
vature contributions depend on the choice of dividing in- 
terfaces. This is actually possible if an additional albeit 
plausible assumption is made (see below). In the next 
step the bicone filled partially with liquid and partially 
with vapor is analyzed. What is found for the curvature 
terms in the previous step for the container filled with 
an homogeneous fluid is then transferred to the new sit- 
uation of a container filled with an inhomogeneous fluid 
without any modification. 
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1. Homogeneously filled bicone 

We decompose the grand canonical potential into vol- 
ume, interfacial, line contributions, and we include curva- 
ture contributi ons due to the curv ature of the solid-fluid 
wall (see, e.g., fl58l [l59L flBcl I 16lj . The curvature contri- 
butions could be also expressed as a curvature correction 
to the solid-fluid interfacial tension. Terms which do not 
scale at least with a linear container dimension are dis- 
regarded: 
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(6.7) 



which is the grand canonical potential for one half of the 
system and p is the pressure of the fluid, vMhe the vol- 
ume of one cone, A^ e _ 8sf the area of the side surface of 

(i) 

one cone, J, is the interface tension of a planar solid- 
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is the 
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gas or solid-liquid interface, respectively, C { 
mean curvature of the side surface integrated over the 
whole area of the side surface of one cone with k , ,•. as 

sg(sl) 

the corresponding thermodynamic coefficient, and r c ^ nt 
is the line (edge) tension associated with the joint be- 
tween the two cones (see Fig. [8]) of length L^ nt . The 
superscript W indicates that geometric quantities as well 
as the interfacial and line tensions and the coefficient of 
the curvature term depend on the choice (i) of the di- 
viding interface between solid and fluid. The integrated 
mean curvature of the side surface of a cone with base 
radius is 
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If we again compare the decompositions for the two 
choices (1) and (2) of the dividing interface shown in 



Fig. [H] and if we use cr , 
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variant of Eq. (|2.3p not restricting the generality of the 
following arguments) and if we neglect terms which do 
not scale at least with rW we obtain the relation 



K sg(sl)°' 



(1) 

cone— ssf 



1 



+ p- 



.(2) 
sg( 



T (l) rW 
cont cont 



,(i) 



K x "' C 
sg(sl) cone — ssf 



-C 1 ) 



■J?) r(l) 
> cont cont 



sin f COS ( 
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-2vrr (1) . (6.9) 



Lcont around the joint of the two cones and attribute to 
this edge a total (integrated) mean curvature based on 
the following arguments. For a container with smooth 
walls without any edge, differential geometry provides 



the general relations (see, e.g., Refs. [162|, LL63() 



y(2) _ y(l) _ M + CS 2 



and 



A {2) = A (l) _ + XS 2 



(6.10) 



(6.11) 



expressing the changes in volume and surface area of a 
container, associated with an infinitesimal parallel shift 
5 of its surface, in terms of the area A, the total mean 
curvature C, and the total Gaussian curvature X. From 
the calculated volume and area changes and from Eqs. 
(|6.10p and (|6. 1 1|) , the total mean curvature for the bicone 
(with rounded edge) turns out to be 



C = 



2irr 



sin 9 cos 9 



(6.12) 



whereas adding up the integrated mean curvatures of the 
side surfaces of the two cones (Eq. (16.8[) 1 gives 



a 



2nr sin 9 



bicone 



sin a cos ( 
This means that a contribution 



2nr cos 2 9 



(6.13) 



(6.14) 



to the total mean curvature of the container is missing. 
Obviously this contribution can be attributed to the line 
L C ont at which the two cones are glued together and the 
curvature attributed to that line can be realized by de- 
forming the surface in its vicinity into one that is diffcr- 
entiable. We use now this observation in order to replace 
the line (edge) term by an equivalent curvature term: 



2 7"cont -^cont y ^ ^sg(sl) Cscam 



irr cos 2 9 

= K sg(sl) — 7 T. 

y ' sin 9 cos 9 

Combining Eq. (|6.9p and Eq. (|6.15|) one finds 



V sg(sl) 



^sg(sl) 



■pS 2 



(6.15) 



(6.16) 



Without a further assumption it is of course not possible 
to obtain the transformation behavior upon shifting the 
dividing interface for t S g(si) and r cont independently. In 
order to proceed we consider a rounding of the edge at 



2. Bicone filled partially with liquid and vapor 

The grand canonical potential of the bicone filled par- 
tially with liquid and partially with vapor (see Fig. [5]) 
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decomposes as 

2 F'cone T ^conc-sg °sg T ^-cone-si °sl 

+ K sg °cone-sg + K s l °conc-sl + ^lg CT lg 

+ ^4o ) „tril + LWrW , (6.17) 

where -A;„„„ scr and A;*'„ Q . are those areas of the cone 
side surface which are in contact with the gas and the 
liquid phase, respectively. <rig and are the interface 
tensions of the planar solid-gas and solid-liquid inter- 
faces. Ccone-sg an d ^cone-si are * ne mean curvatures 
of the cone side surface integrated over that part of the 
surface which is in contact with the gas phase and the 
liquid phase, respectively, nig and K?? are the corre- 
sponding curvature coefficients. A[g is the area of the 
planar liqid-gas interface and eri g the surface tension of a 
planar liquid-gas interface. The last but one term in Eq. 
(|6.17|) . already familiar from the homogeneously filled 
bicone, gives the line contribution from the edge along 
which the two cones are glued together. The last term is 
the contribution from the solid-liquid-gas three-phase- 
contact line. The superscript W indicates which of the 
quantities depends on the choice of the solid-fluid divid- 
ing interface (cr|g — cr? is independent of such choices). 
As in the previous subsection we compare the decompo- 
sitions for the two different choices of dividing interfaces 
indicated in Fig. [SJ we neglect all higher order terms 
which do not scale with a characteristic system size, we 
further use Eq. (|6.9[) which has been obtained for the 
homogeneously filled bicone, and finally end up with the 
relation 

(4 1) ~««) 7 rr 1 «cot0 + 2^rW 

= (^ ) -^) 7 rr«cote + 2 7 rr£V( 2 ). (6.18) 

Equation (|6.18[) states that a certain combination of cur- 
vature and line contributions is invariant under a notional 
shift of the solid-fluid dividing interface. In order to fix 
the transformation behavior of the individual contribu- 
tions under a notional shift, an additional convention is 
required. Insisting on using Eq. (|6.16p . which appears to 
be a most plausible choice for the kind of systems consid- 
ered, fixes the transformation behavior of the difference 
of curvature coefficients which appear in Eq. (|6.18p . This 
leads to Eq. (|3.3[) which tells that r does depend on the 
choice of the solid-fluid dividing interface. 



C. The drop revisited 

The surface tensions at planar solid-fluid interfaces de- 
pend on the position of the dividing interfaces as given 
in Eq. (|2.3p . Therefore, the difference a aj — tr^ 7 of the 
gas-solid and the liquid-solid interface tensions at coex- 
istence of the gas and the liquid phase does not depend 
on the choice of the solid-fluid dividing interface. In our 
analysis of the drop in Sect. V we could make use of that 
property because in defining r we have used solid-fluid 
interface tensions at the very same pressure irrespective 
of whether the fluid phase is the gaseous phase outside 
the drop or the liquid phase in the interior of the drop. 
We have also discussed the possibility of introducing in- 
stead a quantity (<x Q7 — erg-)' in which the two interface 
tensions arc mcassurcd at different pressures. We also 
have given reasons why we have not pursued this possi- 
bility further. However, in view of the conflicting results 
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FIG. 9: A sessile liquid drop on a planar substrate (7). Two 
choices for the substrate-fluid dividing interface separated by 
a distance Sh are shown. The contact angle 9 is only shown 
for that choice of the dividing interface, which is indicated by 
the solid horizontal line. The pressure in the liquid phase (/?) 
differs by Ap from that in the gas phase (a). The drop is a 
spherical cap with radius R — r/ sin 9, area A = 2tvR 2 (1 - 
cos 9), and volume V = (n/3)R 3 (2 - 3 cos (9 + cos 3 9). 

for the transformation behavior under notional shifts of 
dividing interfaces between different definitions of line 
tensions we now seriously consider this option. 
In order to make our reasoning as transparent as possible 
we discuss now only a notional shift of the substrate-fluid 
dividing interface as shown in Fig. [HI but leave the liquid- 
gas dividing interface fixed. We again use the decompo- 
sition in Eq. (|5.9p for two positions of the substrate-fluid 
dividing interface and use at first the relation 

(<r®-<r%)-(<r$-<r$)=0 (6-19) 

which is valid if both interface tensions are measured at 
the same pressure (and if for both the a-7 and the /3—j 
interface the dividing interface is positioned at the same 
height). After neglecting terms which do not scale with 
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at least a linear extension of the system one arrives at 

2nr (r« - r< 2 >) = -Ap (v^ - V^A 
+a aP (R)(AW-AU)+n(a M -a^)((rW) 2 (r«) 2 ). 

(6.20) 

Since in Eq. (|6 . 20[) neither the liquid-gas surface ten- 
sion <J a p{R) ( in the relevant order) nor the difference 
(<7/3 7 — (Ja-y) depend on the choice of the substrate-fluid 
dividing interface, the line tension r is the only physi- 
cal quantity that can pick up all the notional changes of 
volume and interfacial contributions. In particular, the 
notional change of the volume contribution (in leading 
order V^ 2 ' — = —nr 2 Sh and Ap oc 1/r) does also 
contribute to the notional change of r. Evaluating Eq. 
(|6.20p up to leading order, i.e., neglecting corrections to r 
of the order of 1/r , one inevitably comes to the conclusion 
that r is independent of the choice of the substrate-fluid 
dividing interface. 

The only sensible way how an alternative definition of 
the line tension t could be introduced is to use in the de- 
composition scheme of Eq. ([5. 9ft a solid-liquid interfacial 
tension ap 7 (p+Ap) taken at the pressure p+Ap of a bulk 
liquid in contact with the solid, whereas the solid-gas in- 
terfacial tension a aj (p) is taken at the pressure p. One 
may argue that this way one mimics the actual conditions 
at the interfaces of the liquid-drop-vapor-solid system. 
On the other hand the interfacial tensions a ai and crg 7 
are not measurable individually, whereas the difference 
(T/3 7 (p) — a a -y(p) is a measurable quantity. Nevertheless, 
we now seriously consider that definition of t. In order 
to study the transformation of this newly defined r under 
notional changes we proceed as at the beginning of this 
subsection with the only difference that instead of Eq. 
(j6.19[) we have to use now 

(a^(p)-a^(p+Ap))-(a^(p)~a^(p+Ap)) = ApSh. 

(6.21) 

If we carry out the analysis neglecting all terms which 
obviously can only give rise to corrections to r of the 
order 1/r, i.e., terms which cannot contribute to a line 
tension in its strict sense, we arrive at 

r (2) _ r (l) = 

Sh 

— [a af3 (R) + (o$ (p + Ap) - a« (p) ) cos 6 . 

(6.22) 

We now replace af^(p + Ap) in Eq. f!T2"2"|) ) by cr$(p) 
because the difference of the two quantities is of the order 



of 1/r (we then can skip the superscript W in the differ- 
ence between the solid-liquid and the solid-gas surface 
tensions) and with the same argument we replace a a p(R) 
by its planar limit <r^. We then can use Young's law and 
finally arrive at 

r (2) _ T (i) =a ™ fj shsm9 , (6.23) 

i.e., we again find the transformation in Eq. ()3.3j) . as we 
did repeatedly in discussing systems with planar liquid- 
gas interfaces. The reason why this second definition of 
the line tension leads us again to the transformation in 
Eq. (13. 3p can be traced back to the fact that in the second 
definition a contribution to the notional change of the 
volume contribution Ap (V^ — VW) is absorbed by a 
corresponding notional change of a a ~ ( (p) — a 01 (p + Ap) , 
whereas in the first definition this possibility does not 
exist. 

1. Relation between the two line tensions 

We now derive a relation between the two values of r 
determined for the same system under identical condi- 
tions, but using the two different definitions introduced 
above. In order to do so we compare the two decom- 
positions of the grand canonical potential, on which the 
two definitions are based on, for a given, but otherwise 
arbitrary, choice of dividing interfaces. In order to dis- 
tinguish the two definitions we denote with r the value 
of the line tension determined according to the decompo- 
sition scheme introduced in Sect. V (see Eq. (|5.8p with 
a ai (p) and <j/3~,(p)), which we sketched again at the be- 
ginning of this subsection. With r w we denote the line 
tension determined via the second decomposition scheme 
(see Eq. (|5.8p but with a ai {p) and cr^ 7 (p -I- Ap)). On 
the basis of its transformation behavior and because ef- 
fects due to the Laplace pressure, which do not occur for 
planar liquid-gas interfaces, are extracted from its defi- 
nition, we identify r w with the line tension as defined for 
the previously discussed systems with planar liquid-gas 
interfaces (it shall be understood that the identification 
applies to the line tension in its strict sense, i.e., disre- 
garding subleading terms). The comparison leads to the 
relation 

t = r w + T - [a 01 {p + Ap) - a 01 (p)] . (6.24) 

For a one-component fluid this relation can be re- 
expressed by using 

a(p + Ap)=a(p) + ^^Ap , (6.25) 
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where p is the chemical potential, and 

- JL 

where /5b is the bulk density of the fluid and 



p(r)d 3 r - p b V" 



(6.26) 



(6.27) 



whith r as the excess adsorption. Equation (|6.27|) can be 
easily derived from the definition of the interface tension 



a = J (n[p eq (r)]+ P V). 



(6.28) 



where p q(j) is the equilibrium number density mini- 
mizing f2, and the general functional form of the grand 
canonical potential: 

fi [p(r)} = T (p(r)) + J p(r) (F ext - n) d 3 r . (6.29) 

Combining Eqs. (|6.24[) - (|6.29p and expressing the 
Laplace pressure Ap in terms of the liquid-gas surface 
tension one obtains (pb,/3 is the bulk density in the (3 
phase): 



r 



Ph,0 



-a a/3 smf 



r^ 7 is the excess adsorption at the planar interface. 
In Eq. (|6.30p both terms on the right hand side depend 
on the choice of the liquid-solid dividing (f3—j) interface 
but these dependences cancel and thus r is independent 
of that choice as it should. 

The discussion in the present section has shown that the 
contradiction between our statement in Sect. V, that r 
is independent of the choice of dividing interfaces, and 
Eq. (|3.3[) is resolved as follows. We have found that the 
definitions of the line tension at a substrate-fluid-fluid 
interface, which have been chosen on one hand in Sect. 
V and on the other hand in Eq. (|3.3p . are different. The 
line tension, which appears in Eq. (|3.3p , is r w and differs 
from t as defined in Sect. V. The possibility to define two 
different line tensions arises because it is not obvious how 
to subtract the contribution to the grand canonical po- 
tential stemming from the planar /3~7 (liquid-substrate) 
interface of a drop of /3-phase in contact with a substrate. 
The liquid-substrate interface tension ap^ is influenced 
by the Laplace pressure, but <7^ 7 is not accessible. How- 
ever, the difference ap 1 (jp) — tr Q7 (p) of the fluid-substrate 
interface tensions is measurable via the contact angle for 
large drops. If in the decomposition scheme defining the 
line tension we use an interface tension at a pressure de- 
viating from the true pressure in the interior of the drop, 



i.e., (7g 7 (p), in order to avoid the somewhat artificial com- 
bination aj3~ f (p + Ap) — a ai {p), one is led to the definition 
r for the line tension. If instead one uses cr^ 7 (p + Ap) in 
the decomposition one arrives at the definition r w . Cor- 
responding alternative choices do not exist in the case 
of the contact of three fluid phases, since all fluid-fluid 
interface tensions are measurable individually and since 
the dependence of interface tensions on curvature radii 
parametrize their dependence on the Laplace pressure in 
a natural way. By contrast, the geometry of the solid- 
liquid interface is not influenced by the Laplace pressure 
and it stays planar for all values of the Laplace pressure. 

The definition r w for the line tension has the merit that 
interface contributions due to the modification of the f3-j 
interface tension caused by the Laplace pressure, giving 
rise to a contribution proportional to the linear extension 
of the system, are not implicitly included in r w . However, 
in order to use t w in equations determining the contact 
angle, it is necessary to re-express a number of equations 
given for the drop in Sect. V. Equation (|5.44p expressed 
in terms of r w and its notional derivatives turns into 
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dR 



= sin 9<j a r : 



(6.30) an d instead of Eq. (|5.45|) one obtains 



<7 a /3 COS 6 + (apj (p + Ap) - er Q7 (p) ) = — — — sin f 



or in an alternative notation 



(6.31) 



dr w 



dR 
(6.32) 



o- a /3 cos 9+ (af3j(p)-a a <y(p)) = -{a l3l {p+ Ap)-a ^(p)) 

dr™ 



r 



— sm( 



dR 



(6.33) 



The relations between the notional derivatives of r and 
t w are 
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dR 
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2 sin 9 



{<7/3~f(P + Ap) - <T/3<y(p)) 
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2 sin^ 



(c,3 7 (p + Ap) -(jfrip)) 



In terms of the variables 9 and r one obtains 
dr w 



A9 



smf 



'a/3 



and 



a a p cos 9 + (cr/3 7 (p + Ap) - a ai {p)) 



T 



dr w 
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sin t/ cos ( 



de 



(6.34) 

rAp 
(6.35) 

(6.36) 
(6.37) 
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and 



0>y (P + Ap) - cr Q7 (p)) = — — 



dr 



(6.38) 



which replace Eqs. (|5.46|) and (|5.47[> . It is interesting 
to see that in Eq. (|6.38f> the term a a p cos 9 drops out 
due to the identity Eq. (|6.36[) . i.e., this term must now 
be included in the term [^r] ■ This is displayed in the 
relation 



= (Jap COS 9 - -(<T/3 7 (p + Ap) - CT/3j(p)) 



+ 



dr 
dr 



+ 



sm V cos f 



dr 
d9 



(6.39) 



Equations (|6.36|) - (|6.39|) show that using t w instead of r 
spoils the clear hierarchy in the various terms describing 
notional changes. 

For completeness we provide also the transformation laws 
for the notional derivatives of r w : 
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(6.42) 



(6.43) 



We make now again contact with a variational ap- 
proach with the constraint of fixed volume. Of course 
one could stick to the variational approach as discussed 
in Sect. V without any modifications. In that case one 
just might want to rewrite Eqs. (|5.55|) and (|5.56|) in terms 
of t w and its notional derivatives. Alternatively, one may 
slightly modify the variational approach by introducing 
a /3~7 interface tension measured at a pressure p + Ap in 
the same way as in the definition of r w . As in Sect. V we 
further introduce stiffnesses and but we do not 
introduce any further stiffness. (One might contemplate 
to endow the /3~7 interface with a stiffness but there is 
no geometric measure characterizing that interface which 
could be related to such a stiffness; the area is already 
used and the radius of cirumference is better attributed 
to the contact line.) In that way one arrives at 
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One also obtains the relations 
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and the transformation law 
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VII. SUMMARY 



We have discussed some conceptual issues that arise in 
a macroscopic thermodynamic description of the three- 
phase contact of either three fluid phases or of two fluid 
phases meeting an inert solid substrate. We have pointed 
out that the concept of a line tension accompanying the 
contact line has to be used with great care. The concep- 
tual difficulties arise because the interfaces between two 
phases are always diffuse and never sharp. Therefore 
a certain degree of freedom exists in positioning, some- 
where in the transition region, an idealized mathematical 
interface (the so called Gibbs dividing interface) separat- 
ing the two adjacent phases. As a consequence, a similar 
degree of freedom exists in the position of the contact 
line defined as the common intersection of three dividing 
interfaces. 

We have analyzed implications for a consistent descrip- 
tion of the contact line following from the existence of 
these degrees of freedom. For that purpose we have in- 
vestigated two representative systems of three-phase con- 
tact: a liquid lens at a fluid-fluid interface and a liquid 
drop on top of a smooth substrate. Both systems are 
used in experimental attempts to determine a line tension 
t from the dependence of contact angles on the system 
size. We have defined a prescription for decomposing the 
grand canonical free energy of a liquid lens or drop into 
volume, interface, and line contributions which renders a 
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line tension r independent of a particular choice of the 
dividing interfaces. The prescription rests on geometri- 
cal definitions and on the notion that those interfaces 
of a lens or drop system, which are spherical segments, 
should be described in the same way as the interfaces 
of completely spherical drops. In particular this means 
that the pressure drop across the curved interface is re- 
lated to the surface tension (and surface stiffness against 
changes of the radius of curvature R) via the generalized 
Laplace equation (|5.11[) and that the surface tension of 
a spherical interface is independent of a notional shift 
[dR] of the dividing interface up to and including conti- 
butions proportional to 1/R for large R. We have also 
used the fact that interface tensions of planar interfaces 
between two fluids in thermal equilibrium are indepen- 
dent of the position of the dividing interfaces. The same 
is implied for the difference between the substrate-gas 
(substrate-vapor) and the substrate-liquid interface ten- 
sions although the substrate on one hand and the fluid 
phases on the other are not in equilibrium. In fact this 
is true if the substrate-liquid and the substrate-gas in- 
terface tensions used in the decomposition scheme corre- 
spond to the same pressure, a prescription which has the 
advantage of relating the difference of these two quan- 
tities directly to the contact angle of macroscopic drops 
which is a measurable quantity. Our result, that r is in- 
dependent of the choice of dividing surfaces, is first of all 
a generalization to curved interfaces and curved contact 
lines of what is known for the line tension of a straight 
contact line at the intersection of three planar interfaces 
in a genuine three-phase contact 0, Q . A further gener- 
alization is the one from a genuine three-phase contact 
to a contact between two fluid phases and an inert solid 
phase. It should be noted, however, that in the latter 
case an alternative definition of the line tension, denoted 
as t w , is possible which seems to be more natural for sys- 
tems containing planar liquid-gas interfaces and planar 
or curved solid walls and more useful for the purpose of 
computing line tensions by exploiting the simplest possi- 
ble geometries. It should be emphasized that t w does 
depend on the choice of dividing interfaces. The dif- 
ference in the definition of this alternative line tension 
relative to the previous one rests on choosing in the de- 
composition scheme of the grand canonical potential a 
different substrate-liquid interface tension which is not 
taken at the pressure of the gas phase, which is implied 
for the substrate-gas interface tension, but at a pressure 
which is enhanced by the Laplace pressure. We have 
provided a simple relation between the values of the line 
tensions corresponding to the two alternative definitions 
(Eqs. KM and (IQOll L 



We have further pointed out that the generalized Neu- 
mann or Young equations obtained from a minimization 
principle by simply adding a line-tension contribution to 
the free energy suffer from internal inconsistencies. The 
purely geometric relations between two sets of contact 
angles obtained for two different choices of the divid- 
ing interfaces are at variance with those equations and 
with our result from the decomposition of the free en- 
ergy which states that the line tension r is independent 
of such descriptive ambiguities. In the case of a drop on 
a solid substrate using the alternative definition r w of a 
line tension cannot resolve those inconsistencies. 

In order to find equations for the contact angles which 
are internally consistent we have followed two different 
routes. The first one is determined by the observation 
that the grand canonical free energy must be indepen- 
dent of notional (descriptive) changes of the system, i.e., 
of 'parallel' shifts of the dividing interfaces at fixed phys- 
ical configurations. In the mathematical formulation we 
introduced notional changes of the line tension and corre- 
sponding notional derivatives of r with respect to contact 
angle(s) and to the contact-line radius in addition to the 
well established notional derivatives of surface tensions 
with respect to radii of curvature. The second route fol- 
lows a minimization of the grand canonical free energy 
under the constraint of fixed drop or lens volume. In 
contrast to the previous formulations we have included 
into the free energy contributions from the stiffness of the 
interfaces with respect to a change of the radius of curva- 
ture and from stiffnesses of the contact line with respect 
to changes in contact angle(s) and to the contact-line 
radius, all at fixed thermodynamic variables. The equa- 
tions following from these two routes have been compared 
and combined into one (set) of equations. For the lens 
these are Eqs. (|5.38p and (|5.39p together with the re- 
lations in Eqs. ([OT]) - (pT43)) and also Eqs. Ip4| and 
(|5.35p relating the stiffness constants and notional deriva- 
tives. For the drop the corresponding equations are given 
by Eq. (|5.55|) and the relations in Eqs. (|5.56|) and (|5 .51 j) . 
Furthermore we have found that the actual values of the 
stiffness constants depend on the choice of the dividing 
interfaces. The relations between two sets of stiffness 
constants for two different sets of dividing interfaces are 
given by Eqs. (|5.32[) and (|5.33[) for the lens and by Eq. 
(|5.54[) for the drop. (In case of the drop one could also 
use the alternative Eqs. (|Q4|) to (|fT49|) .) 

At this point it is interesting to see that for reasons 
similar to those which compelled us to introduce notional 
derivatives of r with respect to contact angles or stiff- 
nesses of t with re spect to changes o f con tact angles, 
Djikaev and Widom |l39| (see also Ref. |l4fj |) introduced 
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a kind of 'derivative' of r with respect to contact angle 
in order to restore invariance against notional shifts of 
dividing interfaces of their linear adsorption equation for 
a straight contact line at a genuine three-phase contact. 
Eve n clo ser to our discussion is the work by Rusanov et 
al. |l37j |. which differs, however, from ours in a num- 
ber of important points. First, Rusanov et al. discuss 
only drops on a substrate, whereas we discuss both lenses 
and drops. Second, in contrast to Rusanov et al. we 
have included notional shifts of the substrate-fluid inter- 
faces. Third, Rusanov et al. have not used a standard 
variational principle at constant volume (typical exam- 
ples for the use of this principle are given, e.g., in Refs. 

155]) and they did not give transformation laws 
between values of stiffness constants for different choices 
of dividing interfaces. Finally we have tried to make clear 
at which points in our line of arguments there is still 
room for chosing different conventions and what follows 
as a necessity for any sensible convention. 

Next we discuss the consequences of our investigations 
for interpreting experimental data. Before doing so, we 
provide explicit expressions for the change in contact an- 
gles P = 2tt — (a + 7) or 9 with respect to their limiting 
values ao, /3o-> 7o or $0 f° r macroscopicly large lenses or 
drops, respectively (see Fig. 3). From Eqs. (|5.39[) and 
(|5.38p together with Eq. (|5.10[) and the corresponding 
equation for <r^ 7 we obtain 



cos (5 — cos /3 



sin/3o 



25^3 cos a 



+ 25 J 7 cos 70 + 



2 r 



,dr 
dr 



a af3 sin ao + £7^ sin 70 



cos ao dr 



J af3 



+ 



cos 70 dr 







lnr 



(7.1) 



Similarly, Eq. (|5.55[) together with Eq. (|5.10|) yields 



cos 9 — cos 9q = 



rcr 



a/3 



oat „oc dr 



sin #o cos t/ 



dr 
dr 







lnr 



(7.2) 



The leading type, =pf j of correction terms arises due to 
algebraically de cayin g dispersion forces among the par- 
ticles (see Ref. ^M)- In E q s - CEO) and JE21 r is the 



radius of the circular contact line, cr|£ the interface ten- 
sion of the planar £—v interface, 6j v the Tolman length of 
the £-1/ interface, ^\ etc. are stiffnesses against changes 
of contact angles or of the radius of the line, which are 
attributed to the line, and r is the line tension. In the 



traditional analyses of size dependent contact angles only 
the term proportional to t is included. However, the 
omitted terms give contributions to the right hand sides 
of Eqs. (|7.ip and (|7.2p which are comparable in magni- 
tude to that term. Although r is independent of the cho- 
sen dividing interfaces, the stiffnesses ^| and in the 
case of the lens and ^| and ^ in the case of the drop 
are not. Therefore, the additional terms in Eqs. (|7.1|) 
and (|7.2p containing these stiffnesses even depend on the 
positions of dividing interfaces which may be chosen ar- 
bitrarily within a certain range. This is in accordance 
with the fact that this dependence is also present for the 
contact angles on the lhs of Eqs. (|7.ip and (17. 2p (see step 
1 in the protocol given below). Moreover, the changes in 
the values of these particular terms with the position of 
dividing interfaces are as big as the term proportional to 
r itself. In order to demonstrate this we compare r with 
7T- 1 , :n- 1 , t£ I , or ^ I as well as with the terms of the form 

da 1 ' d7 1 ' d9 1 ' dr 

25 T (J involving the Tolman length S T . The change, e.g., 
of jp| with a shift of the a-(3 dividing surface by [di?i] 
is equal to — <7 Q/ 3[di?i] For typical values of a a/ 3 of the 
order of 10 _2 J/m 2 , and for [di?i] of the order of lnm 
one obtains a change in the value of 4r | of the order of 
10~ n J/m, which is just the typical value of calculated 
or 'measured' line tensions. Similar estimates hold for 
the other terms in particular also for those of the form 
25 T <J if the reasonable assumption is made that 26 T is 
not much smaller than 1 nm (typical values of 25 T ob- 
tained theoretically are o f the order of half a molecular 

diameter, see, e.g., Refs. [HEEii, [Ti3i)- 

In order to analyze experimental data with the help of 
Eqs. (|7.ip and (|7.2p a certain protocol has to be followed 
as discussed below: 

1. Giving experimental values of contact angles (3 or 
9 requires to define the dividing interfaces relative 
to which the contact angles are measured. (So far 
this information is missing for basically all pub- 
lished experimental efforts to measure the line ten- 
sion.) After the dividing interfaces have been cho- 
sen (iso-density surfaces) , the spherical parts of the 
interface profiles of the drop or lens have to be de- 
termined by fitting spheres to their central parts 
(e.g., to the top of the drop). In case of the drop, 
one such sphere intersects with the planar solid- 
fluid interface (chosen by convention) at the three- 
phase-contact line. The contact angle at this line is 
obtained by a tangential-plane construction to the 
sphere, i.e., to the spherical extrapolation of the 
liquid-gas interface down to the solid-fluid inter- 
face. In the case of the lens, two spheres are de- 
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fined by the central part of the lens. These spheres 
intersect at the three-phase-contact line, at which 
tangential planes to the spheres define the contact 
angles. 

2. The dependence of these angles on the lens or drop 
sizes, as characterized by r, has to be studied in 
the leading asymptotic behavior for large r. This 
dependence is not only determined by the line ten- 
sion t calculated for a straight contact line but it 
depends on additional material parameters. There- 
fore, the common practice of infering t from f3(r) 
or 9{r) according to Eqs. (|4.3|) or (|4.4|) is not valid. 

3. If one would compare the rhs of Eqs. (|7.ip and 
(|7.2p with theoretical quantities, these would have 
to be evaluated for that choice of the dividing in- 
terfaces for which the experimental data on the lhs 
are given. 

4. The values of certain stiffness constants appearing 
on the rhs of Eqs. (|7.ip and (|7.2[) depend sensitively 
on the choice of dividing interfaces. However, the 
transformation behavior of these material param- 
eters between different such choices is known and 
given by Eqs. (|5.32p and (|5.33p for the lens and Eq. 
(|5.54j) for the drop. 



All quantities on the rhs of Eqs. (|7.ip and (17. 2p are 
accessible to theoretical computations, e.g., based 
on density functional theory. Methods that can be 
used to calculate the Tolman le ngths are described 
in the literature (see, e.g., Re fs . 
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liifiiifiii. 



In calculations of the line tension much care has to 
be taken not to pick up additional artificial line or 
edge contributions. In the likely case that within 
a theoretical approach one has calculated r w , its 
relation to t is given by Eqs. (|6.24p and (|6.30D . 
provided both quantities apply to the same thermo- 
dynamic conditions. 

Once the density profiles for a lens or drop have 
been computed, the stiffness constants etc. 
may be determined from their relations to the no- 
tional derivatives of t (see, e.g., Eq. (|5.56p ) or in 
the case of 4j | also to the stiffness constant of the 
interface (see, e.g., Eq. (I5.5ip ) which again is given 
by a notional derivative (Eq. (|5.13p ). Similar equa- 
tions which apply to the case of the lens are Eqs. 
(pQgjl . (IQ51) . rfOTjl . {Q2jl . and (J5~43]). r itself as 
well as the notional derivatives of r follow from de- 
composing the grand canonical potential for a series 
of different choices for the dividing interfaces. For 



instance one might decompose the grand canonical 
potential fl of a given lens for a series of different 
values of [Ri] and [Rz]. Carrying out this decom- 
position up to the order 1 jr and making use of pre- 
vious knowledge about the notional derivatives of 
the surface tensions (see Eq. f|5.12[) ). the notional 



derivatives 



At 



and 



dr 



dR 2 



can be determined 



from the dependence of r on [Ri] and [R2]. These 
notional derivatives may be then converted into the 
notional derivatives [4^1 and 



d/V 2 



sin a cos 



(g+Tp 



r sin(a+7) 

sin a [" dr 1 
r sin(a-f 7) L da J 



dr 1 
.da J 



r sin(a+7) 
sin 7 cos(q-|-7) 
r sin(a+7) 



via the relations 
and 



larly the notional derivatives and can be 
determined in case of the drop and converted into 
[%] and [g] via the relations [&] = ^ [%] 

[&] and [&] = 



"drl 



cos 8 [" dr 1 
sine Ldhl 



Although each individual quantity entering the rhs of 
Eqs. (|7.ip and (17. 2p can be calculated separately in the 
way indicated above, measurements of the size depen- 
dent contact angles of drops or lenses can provide only 
certain combinations of material parameters. Their val- 
ues depend sensitively but in a known way on the choice 
of dividing interfaces. 

There still remain a number of open questions. 

1. How can the various stiffness constants attributed 
to the contact line be measured? The answer to 
this question requires the extension of the princi- 
ples developed here for simple geometries to more 
general geometries. How do the stiffness constants 
of the contact line influence the equilibrium shapes 
of more complex structures, e.g., of liquid bridges 
between solid substrates or of lenses and drops dis- 
torted from their ideal shape due to external forces? 

2. How can the stiffness constants of the contact line 
introduced here be related to the 'derivatives' of 
r with respect to contact angl es wh ich were intro- 
duced by Djikaev and Widom 1391 ] in order to re- 
store invariance of their linear adsorption equation 
against notional shifts of dividing interfaces? 
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